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ABSTRACT 


This study provides information on the distribution 
of surficial sediments, and morphological characteristics 
of the Elbow River. basin, in order to facilitate future 
research in this area. Discussion of the evolution and 
present behaviour of the Elbow River, and notable events 
concerning the development of the landscape is included. 

Sediment samples and morphological information were 
collected in the field. The results of the laboratory 
analysis of grain size were analysed statistically using 
principal component analysis. The morphological information 
was supplemented by mapping from aerial photographs. Results 
included a morphological map and a map of the surficial 
geology. 

The conclusions of this study are: 

1. Principal component analysis can be used to 
discriminate between sediment types. 

2. Neither antecedence nor superimpositions are 
conclusively proven to be responsible for discordance between 
the Elbow River and geological structure in the foothills 
region of the study area. 

3. The distribution of the surficial material is 
closely related to limits of advance of Cordilleran and 
Laurentide ice masses. 


4. The present Elbow River. is in a state of 
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quasi-equilibrium, due to a decrease in mean annual flow 
Since neo-glacial times which has resulted in a change in 
the sediment transport characteristics. The result is 
adjustment of slope in the longitudinal profile to; (a) 
changes in the amount and size of sediment transported; and 
(b) the lowering and/or removal of the influence of 
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INTRODUCTION 


This study concerns aspects of the geomorphology of 
the Elbow River basin (Fig.l). In particular, the work 
studies the basin's morphological and sediment characteris- 
tics, and analyzes the behaviour. of the present Elbow River. 
The scope of the. study is limited in extent and nature by 
time available for both fieldwork and laboratory analysis. 
The major objective is to provide a basis for assimilation. 
of earlier, more specific studies and to indicate useful 
future research. 

At this time most of the studies relating to the Elbow 
River and its basin fall into one of three groups. The first 
group contains geological studies largely concerned with oil 
or coal exploitation. The second group of studies involves 
specific aspects of the hydrology of the basin. The third 
group of work relates chiefly to the stratigraphy of the area. 
This latter group usually includes the Elbow basin only in 
terms of studies centred on the geomorphology of the whole 
Bow River valley. Specific reference to each of these groupe 
of studies is made below. 

The earliest references to the Elbow River are nine- 
teenth century evaluations of the natural history of the area 
(Dawson, 1885). In the early twentieth century, as mineral 
exploitation increased, more geologically oriented reports 


began to appear such as Cairnes (1907), Hume (1928), Slipper 
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(1921). By 1930 more detailed geological studies were carried 
out by Hume (1930 and 1931). The most complete geological 
inventory of the area was prepared by Beach (1943) and pub- 
lished by the Geological Survey of Canada. 

In 1933 the Glenmore Reservoir (Fig. 2), a dam on the 
Elbow River which still supplies water and hydro-electric 
power to the City of Calgary, was constructed. Public aware- 
ness of the Elbow River and its behaviour increased; previously 
only a few individuals who farmed the. low, rich, flat land 
upstream of the reservoir were directly concerned with the 
hydrology of the river. It is of interest to note that the 
name of this river was taken from the description the early 
farmers gave to the channel pattern (Pers. comm., Mr. Murray, 
Calgary). The erratic behaviour of the river has caused 
several inquiries, involving cases of individuals who built 
or farmed on the floodplain and incurred loss or damage of 
property due to channel changes or flooding (Pers. comm., Mr. 
A. B. Hollingshead, Research Council of Alberta). 

The presence of the reservoir created a need for 
reliable information concerning flow duration and the source 
of water. Hence, discharge recording stations have been in 
operation for over forty years in the basin. In 1961 Meyboom 
carried out a study on groundwater resources for the City of 
Calgary. The Elbow River had low groundwater seepage rates 
and consequently a high runoff fraction of total precipita- 
tion. In 1965, Deeprose reported on the Elbow River floods. 


There are two studies of sediment infill rates for the Glenmore 
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Reservoir by Yaremko and Mustapha (1967), and Hollingshead 
(1969). Concern regarding the behaviour of the Elbow River 
and the lack of knowledge of sediment transport rates for the 
river upstream of the reservoir, led to a study by Hollingshead 
(1967 and 1968). This study measured bedload as well as sus- 
pended load for selected reaches downstream of Bragg Creek 
discharge station (Fig. 2). This work by Hollingshead (1967) 
required the building of a third discharge station at Beaver 
Flats (Fig. 2). 

At present the Elbow River basin: is used.extensively 
for recreation by Calgarians. This pressure on the area has: 


caused concern about pollution of Calgary's water. supply. 


Because of future possible coal mining activity in the area 
the Eastern Rockies Forest Conservation Board agreed to 
undertake a study of the capacity of the Elbow River to carry 
potential waste products from strip mines. The study was 
carried out using artificial tracer materials and determined 
the optimal waste-bearing level that the river could sustain 


under varying flow conditions. A pilot study had indicated 
the significance of dispersal of waste as.the rising limb of 
the peak flow curves were developing. (Pers. comm., Mr. W, 


Poliguin. of the E.R.F.C.B.). Finally, the forest.rangers are 
becoming more aware of despoliation of the forestry area 
caused by modern types of vehicles which have access to 
previously inaccessible areas and by commercial interests 


which cause land use conflicts. 


The topics studied by others in the Elbow River basin 


have been highly selective. There exist few comprehensive 
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studies of the physical environment other than of factors 
pertaining to the forestry environment. At the same time, 
as shown above, there appears a need for. more. general research 
in the area and especially for research which can. be used in 
the compilation of an inventory of the basic factors, of the 
physical environment. 

The detailed studies which refer to the study area 
are mentioned in the following chapters. The lack of research 
into the morphological and sediment content of. the basin is 
evident in the general manner in which the Elbow River basin 


is referred to in much of the literature. reviewed. 
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CHAPTER I 
THE STUDY AREA: A DESCRIPTION 


The Elbow River is a tributary of the Bow River, with 

the confluence in.the City of Calgary, Alberta. (Fig. 1). 

The headwaters of the Elbow are in the Rocky Mountains, 

seventy miles to the southwest of Calgary. The river traverses 
the eastern ranges of the Rocky Mountains, through both the 
Alberta foothills and the western edge of the Prairies. 

The geological boundary between the Rocky Mountains 
and the foothills is usually taken as the major thrust sepa- 
rating Palaeozoic from Mesozoic rocks (Beach, 1943). Within 
the Elbow basin the thrust is found on the eastern flank of 
Nihahi and Compression Ridges which run NW to SE (Fig. 20). 
Topographically, the foothills are found to the east of a 
line running NW to SE through Beaver Flats. This line sepa- 
rates the rugged, glacially denuded landscape to the west 
from a region of increasing glacial infill to the east. 
Relief is lower to the east and slopes are less rugged. The 
tree line is found at an approximate elevation of 7500 feet. 
As is typical of most of the Rocky Mountains, soils are mainly 
azonal. In the foothills the soils are chernozems, and the 
ree relief gives to a heavy forest cover with lodgepole 
pine grading to the east into poplar. The writer estimates 
the tree cover to be at least 85 to 90 percent for the foot- 


hills portion of the area studied. A few miles east of Bragg 
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Creek the relief sharply decreases. The lower Elbow basin, 
east of Last Break Hill (Fig. 2) is a broad grassland area 
and soils are more characteristically dark black earth in 
type (Wyatt, Newton and Bowser, 1960). The relief and 
geology of the whole area are discussed more fully in Chapter 
Three. 

The specific study area (Fig. 2) was limited to the 
241 square miles (area taken from basin catchment figures 
given in Hollingshead, 1968) which is downstream of- Beaver 
Flats, It was felt that this size wouid give an optimal 
solution to the two problems of having a large enough area 
to produce meaningful results at the scale desired and having 
a fixed amount of time to travel over that area, collect 
samples, and carry: out laboratory analysis. 

The methods used in carrying out the study of the 
sediments and the morphology mentioned in the introduction 
will now be discussed. Based upon. an initial evaluation of 
the area from aerial photographs maps, and from the study 
carried out by Hollingshead (1967), two important. factors 
could readily be ee ets The first was the apparent sig- 
nificance of geological control first suggested by Beach 
(1943). The second factor was the significance of the content. 
of the surficial deposits marginal to the main valley of. the 
Elbow in the assessment of sediments in, or deposited by, the 
present river. This idea was first suggested to the writer 
by Mr. Hollingshead (Pers. Comm,, 1970). By studying these 


factors it should be possible to determine adequately the 
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genesis, morphology, and present day dynamics of sediment 
distribution in the valley of the Elbow River. 

Most of the time devoted to fieldwork was. spent 
collecting samples of the surficial deposits, although 
observation of morphological characteristics of the landscape 
was made throughout the area. The latter information was 
used in the derivation of the major significant morphological 
interrelations postulated for the area. The compilation of 
this data necessarily requires study and appreciation of the 
Pleistocene and post-Pleistocene chronology of events in the 
area and its hinterland. No attempt is made to establish any 
absolute chronology for the area, though it is necessary for 
clarity to place events relative one to another in time. 

When the surficial sediments are mapped on the basis of 
samples analysed, it is then possible to interpolate and 
consider distributions of sediments over the whole area. This 
enables morphological analysis to be carried out and thereby 
finally yields clues to the river's behaviour in the past. 
This behaviour will primarily relate to general channel 
patterns and to estimates of valley base material and their 
relationship to the morphology and sediments found. 

In carrying out the field and laboratory work, several 
unexpected events were encountered and it may be of benefit 
to future workers planning research in this or similar basins 
to itemize these. A flood of unusual magnitude washed out 
the bridge near the Elbow Ranger Station and this necessitated 


collecting more samples downstream of Beaver Flats than had 
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10 
been originally planned. The existence of the Sarcee Indian 
Reserve imposed rigid restriction on sample collection in the 
centre of the study area. Altimeter traverses were considered 
unreliable because of the rapid daily pressure fluctuations 
in this foothill region. A final restriction was that the 
time required to complete laboratory analyses was considerably 
underestimated, and this aspect of the work was. therefore 
limited to grain size analysis, 

Despite these difficulties nearly all of the objec- 
tives of the study were met, and although the validity of 
the conclusions would have been greater in an area without 
these problems, it is felt that the contributions made by 


this work are nonetheless worthwhile. 
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CHAPTER II 
GRANULOMETRIC ANALYSIS 


This chapter is concerned with sediment sampling, 
mechanical analysis of samples, and the analysis of the grain 


Size data. 
Collection of Samples 


Although random sample plans were prepared, the size 
of the area and the heavy forest cover (Plate 1) forced the 
rejection of this method. Since sites which could be visited 
were Spatially limited, a pragmatic sampling design was used 
and samples collected wherever possible. It is felt, there- 
fore, that the number of samples of each sediment type are 
biased. The number of samples of each sediment type is shown 
in. Table 1. ° Ingaddition, thewirregulagyand oftengexcessive 
distance between sample points (see Fig. 2), combined with 
the scale of variation discovered, in some areas preclude 
reliable interpolation. 

The final sample pattern (Fig. 2) reflects the access 
routes in the basin and emphasizes the main valley base area. 
Notable unsampled areas are the Sarcee Indian Reserve and the 
inaccessible area upstream of Beaver Flats. 

At each site samples were taken as far below the sur- 
face as possible to avoid suspected periglacial interference 


and to minimuze the effects of weathering. Periglacial action 
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TABLE 1 


NUMBER AND TYPES OF SEDIMENTS 


COLLECTED: IN. THE: FIELD 


Total number of sites 


Samples collected (Based on subjective 
field allocation) 


Cordilleran. Till 

Laurentide Till 

Coarse Fluvial Gravels 

Deltaic Deposits 

Fine Grained Fluvial Sediments 
Ash 


Lacustrine 


Total 


a) 


52 


268 
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14 
is evident in the form of ice wedge pseudomorphs and con- 
tortions of sediments in which the structure is apparent from 
the homogeneity of the grain size contents and colour within 
each contortion. Sampling from such structures would not be 
representative of the sediment type. A further precaution 
included the omitting of exposures at road cut locations where 
the road side of the cut was banked (fig. 3). This was because 
of the common practice of highway engineers of ploughing over- 
burden into the cut from the cut side to produce the bank. 

At sites where suites of different sediments were 
found, samples were taken of each of the clearly discernable 


types. 
The size of the sample taken from each sediment type 


at each site was in the range 500 to 1000 grams. In the case 
of sediments thought to be till, the sample was taken in such 
a way as to avoid pebbles greater than 32mm. This was because 


the size of sample required to generate statistically feasible 
results for ReCeEn Ty size content precluded such an under- 
taking in this study. The coarse fluvial samples in which 
boulders greater than 256 or even 512mm in diameter occurred, 
regularly necessitated larger samples being taken and the 
use of a different procedure for collection in- the field. 

The samples taken for laboratory analysis were in the range 
5000 to 10,000 grams and sieved for material less than 32mm. 
The coarse fraction was measured from photographs which were 
taken of each coarse fluvial samples (Plate 2) using a 
plastic one-foot ruler for scale. A correction factor could 


have been applied had it become necessary to distinguish 
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between different types of coarse gravel. The sites (numbers 
83 and 3, shown in figure 2) for which the analysis was 
completed were chosen such that the procedure could be applied 
to both the horizontal and vertical sections (Fig. 4.). 

At each site the top soil was cleared away completely. 
Then a hole was dug to a depth of twice the length of the A 
axis of the largest pebble. The hole was square and the 
length of each side was four times the A axis. This procedure 
was devised by Kellerhals and Bray (1970) to ensure that the 
largest pebble contributed less than three percent of the 
total sample weight. The material removed was sieved in the 
field through 8, 16, 32, 64 and 128mm sieves. The base of the 
hole was levelled and cleared with a soft brush. A grid 
with 100 intersections each 50mm apart was placed over the 
cleared surface and photographed (Plate 3). Finally, the 
pebbles from underneath each intersection were removed and 
the A, B and C axes measured. Material less than 8mm was. 
designated fines. The problem of more than two intersections 
on the same pebble was solved by re-locating using randomly 
selected coordinates. This was done only once for each such 
situation. This procedure was the same in the horizontal and 
vertical case except that in the former, gravels less the 1/2" 
deep were not counted and the underlying pebble measured 
instead. This was done because of the difficulty of arriving 


at an undisturbed and clear surface on the horizontal. 


Mechanical Analysis of Grain Size 


The laboratory analysis was preceded by a grouping 
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19 
of the samples to facilitate a more efficient granulometric 
study. The ash samples, initially determined by their 
surface location at each site and their distinctive orange 
or yellow colour, were analysed mineralogically under the 
microscope with the assistance of Dr. J. Westgate of the 
Geology Department at the University of Alberta. The conclu- 
sion was that the ash samples were clearly identifiable and 
need not be further analysed (Westgate and Dreimanis, 1967; 
Smith and Westgate, 1969). 

The remaining seven groups were sieved. The procedure 
used in the laboratory is given in the flow diagrams of figures 
5 for the fines, and 6 for the coarse fluvial and deltaic 
samples. The two procedures were used since the coarse fluvial 
and deltaic sediments were already clearly distinguishable 
on the basis of pebble size, bedding structures, pebble shape 
and morphological form. The procedure for all other samples 
was designed to allow statistical comparability between each 
group. Hence, sieving was cut off at the upper and lower 
limits of 8mm and 0.063mm respectively. The presence of 82 
samples of at least two till types created the need for the 
upper limit. The distinctiveness of the coarser than 8mm 
fraction could only be truly ascertained by field sieving of 
very large samples. The lower limit was imposed because of 
a shortage of time for the pipette analysis of so many samples. 

The laboratory sieving of samples was carried out at 


half phi intervals. 
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FLOW DIAGRAM FOR LABORATORY 
ANALYSIS OF GRAIN SIZE 
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Analysis of Grain Size Data 


This section is considered as an experiment and is 
presented under the headings of theory, method, results and 


conclusions. 


Theory 

The identification of a sediment sample in terms of 
its mode of genesis is frequently based on several measured 
properties of the sample; for instance grain size, petro- 
graphic and mineralogic content, pebble orientation and 
locational data. At present there is an increasing awareness 
of a need for information concerning the spatial and temporal 
distribution of any one of these properties of a sediment. The 
problem of characterizing a sediment in terms of genesis and 
distribution are complex when only one of the properties are 
considered. This is true of grain size classes, since this 
is. the basic component of all clastic sediments. In this 
Study grain size is the primary property of each sediment 
sample collected. This property is used to characterize 
sediment types. 

In granulometric analysis there is an increasing 
amount of information concerning the characteristics of grain 


size distributions for sediments of various genetic origins. 


Krumbein (1938) has pointed out that grain size analysis is 
largely a problem of statistics. It is not surprising there- 
fore, that most granulometric analyses presently undertaken 
by geologists and geomorphologists are based on the most 


reliable statistical analysis evolved primarily for such work. 
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These techniques are referred to as (1) moment measures and 
(2) graphic statistics. These statistics have been widely 
discussed by such workers as Sahu (1964), Link (1966), 
Shepard and Young (1961), Greenwood (1969), Freidman (BIC L-; 
1962), Landim and Frakes (1968), Folk and Ward (1957), 
McCammon (1962) and Inman (1952). Folk and Ward (1957) have 
developed graphic statistics for grain size analysis to a 
point where most workers now recognize the technique as 
reliable as moment measures (McCammon, 1962; and Middleton, 
1962). The graphic statistics have the advantages of ease 
of computation and their visual graphic properties. 

The graphic statistics of Folk and Ward (1957) are 
based on a distribution approximating a logarithmic- 
probability distribution (Krumbein, 1937; Freidman, 1962; 
Middleton, 1962). The most characteristic divergence from 
this relationship, has given rise to what is considered a 
basic phenomenon of the grain size distribution, namely the 
curve away from the straight: line at’ each extreme ofthe 
distribution (Tanner, 1958). Because of this, one prerequi- 
site of this analysis is that Less” than” five™percenc=or “tne 
individual sample weight is included in the smallest and 
largest sieve class. This graphic analysis which generates mean 
standard deviation, skewness and kurtosis, is based on percen- 
tiles. This means that knowledge of the complete range of grain 
size classes in the sediment is necessary. Therefore, when there 
is doubt concerning the whole range of grain size classes (an 


incomplete data record) or the smallest measured sieve size 
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represents greater than five percent of the individual: sample 
weight, this form of analysis is not justified. In this case- 
another form of analysis is required. 

The descriptive elements generated by graphic statis- 
tics are presently being used to identify and discriminate 
between sediment types (Freidman, 1961; Shepard and Young, 
1961). The basis for this work was outlined by Folk and 
Ward (1957) who wrote: 

"Once a relationship is established in an ideal case, 
where the changes are laid out before the observer in 
their most perfect form, one soon learns to recognise 
the same relationships in less ideal examples, where the 
changes are obscure. The obscure examples, hitherto 
unfathomable, are explained in the light provided by the 
ideal examples." 

This suggests a two phase operation.. The first phase 
is to deduce certain explanatory facts ahi cnmieeinetien 
between sediment types in the ideal case. The second phase 
is the induction from these explanatory facts for the obscure 
case in order to produce the saneudiseraninettene Hence, 
when the four elements (or dimensions) of graphic statistics 
are available, this approach to sediment discrimination may 
be used for samples from a particular area (the obscure. case), 
and tested against the known results for the same analysis 
published by other workers (the ideal case). The discrimina-— 
tion considered above is based on uniformitarianism concern- 
ing genesis of a sediment. The assumption is that one sedi- 
ment type is, more or less, similar in grain size content to 


another. Yet grain size is but one variable which may be 


used to partially describe a sediment. Therefore, an implicit 
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assumption of the above is that all variables other than grain 
size act in the same way for a given sediment type, and that 
the degree of dependence of grain size on all:other variables 
is the same at all times. Greenwood (1969) pointed out that: 
"If a meaningful uniformitarianism can be accepted any 
likely errors focus, at least in part, on the complexity 
of the situation and the great number of variables 
involved. The unrestricted use of any single variable 
as a criterion for environment discrimination is often 
suspect because of its possible interaction with a number 
of other variables. The problem of discriminating between 
sedimentary environments on the basis of measurable 
sediment parameters resolves itself, therefore, into an 
application of multivariate statistics." 
The most unique property of the grain size variable 
is that it may be used to describe any sediment. It is to be 
expected, therefore, that grain size content is commonly used 
to reflect the behaviour of other groups of variables than 
any other single variable. The unknown factor is the degree 
of influence of the 'other variables', The implicit assump- 
tion of Greenwood (1969) is that variance in the grain size 
distribution when measured by grain size classes will be 
manifest in the differences between each grain size class for 
sediments of different genesis. By analysing these differences 
using the multivariate statistic of linear discriminate 
analysis, Greenwood (1969) established that: 
"Multivariate discriminate analysis procedures can provide 
a predictive model for environment discrimination using 


sediment parameters. Samples of unknown origin can be 
classified with a defined probability of. misclassifica- 


tion. 


Problems Encountered in the Analysis 


In this study a large number of sediment samples were 
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collected from a variety of sediment types. The first problem 
was to verify the original diagnosis of sediment type allocated 
to each sample collected on the basis of contents, morphologi- 
cal form and location. The second problem was created by the 
limited time available for laboratory analysis, and involved 
finding an alternative statistical procedure for analysing and 
discriminating between samples of different sediment types. 

The incomplete range of data for each grain size 
curve precluded the use of graphic statistics and moment 
measures. However, Since a constant range of grain size 
classes was used in the mechanical analysis for each sample, 
a basis for the use of a multivariate statistic had been 
provided. Greenwood (1969) applied linear discriminate 
analysis to the four moment measures frequently used to 
describe the grain size curve. Greenwood (1969) used these 
four variables as a means of reducing the original number of 
variables (the grain size classes). However, it was pointed 
out by the same writer that principal component and factor 
analysis could be used for the same purpose. Greenwood 
pointed out that the problem of factor analysis is in the 
interpretation of the new variables which may well distinguish 
between the different sedimentary environments, but prove dif- 
ficult to interpret in terms of their geological significance. 

Within the vast field of linear algebra and group 
theory, Gould (1969) has illustrated the potential use of 
eigenvalues and eigenvectors in the multivariate classifica- 


tory problem. Gould (1969) also points out the usefulness of 
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principal components analysis in the context of the matrix 
properties of eigenvalues. In addition, he points out that 
the parameter being used may have a muiti-dimentional char- 
acter so that the classification cannot be established on one 
scale. The taxonomic properties of principal components were 
illustrated by Kendall (1939) for crop productivity in England. 
The problem of spatial uniformity and/or continuity of a 
parameter is much reduced when the taxonomic theory applied 
to a parameter yields a single scale upon which the classifi- 
cation is based. Gould (1969) suggests that principal compo- 
nents is far more powerful than ratio, interval, binary or 
ordinal scales, and may therefore be considered potentially 
far more useful when a multi-dimentional, spatial, temporal 
and multivariate classificatory problem is encountered. 

In this study the major problem is classificatory; 
the data multi-dimentional and sensitive to multivariate 
dependences, and dependent upon spatial and temporal aspects. 
It was felt therefore that principal component analysis was. 
the most appropriate statisitical procedure to use in order 
to verify the original sediment sample identification. 

The author is unaware of a previous application of 
this technique to grain size data. Therefore, the analysis 
using principal components is presented as an experiment, in 
which the behaviour of each sample and the clusters of samples 
in each component of a particular sediment type are 


investigated. 
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Method 

The analysis was applied to Cordilleran till, Lauren- 
tide till, lacustrine, fluvial sands, and ice wedge material, 
all samples of which had been sieved within the same range of 
0.063mm to 8mm. 

The calculation of percentage and cumulative percentage 
weights as well as the plotting of this data was programmed 
for the computer by the author (see appendix 1 for program 
listing). The plots of the data were made in phi units to a 
linear scale (a semi-logarithmic plot). In addition, the 
cumulative percentage curve was plotted a second time with 
phi units against the logarithm of the cumulative percentage 
(2 x logarithm plot). Figure 7 shows samples of the plots... 

To differentiate quantitatively between the five 
groups consider the data as a matrix (MO-1) of rows of obser- 
vations, and columns of the percentage contribution of each of 
fifteen grain size classes to the sample population. Then 
the method is-to test the significance of.the percentages less 
than 0.063mm for each sample against the original diagnosis 
of sediment type. Secondly, the latent roots, or eigenvalues 
of the matrix (MO-1) are generated. Finally, the component. 
listings of the factor loadings are checked for clusters of 
the five original sediment types, and the existence of any 


new sediment groups. 


Step 1 
The first step is to simplify the matrix (MO-1) so 


that description of each sample is made easier. This is 


as 
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30 
achieved by reducing the original number of grain size classes. 
(15) to seven new classes which approximate the Wentworth 
(1922) scale. The following are the classes used. 

< 0.063mm Clay-silt 
0.063mm to 0.125mm Fine sand 
0.125mm to 0.354mm Medium sand 
0.354mm to 1.000mm Coarse sand 
1.000mm to 2.830mm Gravel 
2.830mm to 4,000mm Gravel and Granules: 
>4.000mm Granules 
Before investigating the new matrix (MO-2) of: grain 
size classes and observations, this matrix is 'rank ordered' 
according to the clay-silt fraction. A computer program was 
prepared for this purpose (see appendix 2 for the program 
listing). The results of the rank ordering of the matrix 
(MO-2) are such that each of the original rows of observations 
‘in the matrix is arranged in an order so that the first column 
(of clay-silt percentages) jis ranked from highest to lowest 
value in the resultant rank order matrix (MO-3). This matrix 
is presented in appendix 3. 


At this point the matrix (MO-3) is studied, and the 


contents of the matrix described. 


Step 2 


The original data matrix (MO-1) is used in this step 
and the multivariate statistics of principal component analysis 
is applied. 

According to Anderson (1958) the eigenvalues and the 
components may be generated from a covariance matrix when the 
original data are in a standardized form. In this study, 


because the data matrix (MO-1) consists of ratio numbers 
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(percentages) obtained from the same dimensional units, this 
matrix may be considered standardized. This procedure is 
useful since a convenient test of significance of the results 
is available in.-this case. However, for convenience, a 
package computer program was used (Dixon, 1968). As in most 
cases, programs which are prepared for this analysis automati- 
cally centre and normalize the original data, and the results 
are then generated from the correlation matrix. This is a 
more efficient computational procedure, but has the disadvan- 
tage that significance testing of the results in much less 
reliable and far more complicated than for the procedure using 
the covariance matrix. 

The computer program (Dixon, 1968) is based on the 
generalised equation, 

(Ce= RE) b<=50 
where b is the matrix whose columns are 

the corresponding eigenvectors to the scalar or eigenvalue }. 
I is the identity matrix. Subsequently the program generates 
[Zz] the vector of standardized components of the model, 


[Zz] = ([W]) . [y'] 
where. 


[W] is the orthogonal matrix of coefficients and [y'] is the 
vector of standardized components given. by 
[y"] [yi~ 
where [y] is given 
by the solution of, 
oes 


where tule is the matrix of covariance. 
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terms (King L.J., 1969). Finally, the program rank ordered 
the column vectors [Z]. Appendix 4 provides a more detailed 
explanation of principal component analysis. 

The results generated by the principal component 
analysis are investigated in the following manner. 

The correlation matrix is checked for groups of 
variables and the degree of interaction of- each variable 
one to another. 

The eigenvalues and percentage of total variance 
explained by each component are checked in order to determine 
which components are most likely to produce meaningful and 


interpretative results. The arbitrary.criteria.for this step 


are that the eigenvalues be greater than 1.0 (King, 1969), 
and the level. of explained variance contributed by each 
component be greater than or equal to 10 percent. A. component 
Satisfying these conditions is. considered acceptable for 
interpretative purposes. 

The third step involves looking for groups: and 
combinations (clusters) of samples within the component.1lists 
which are acceptable for interpretative purposes in the 
previous step. The method of. checking the component.lists 
involved the following procedure. 

The columns of. factor loadings in each component were 
arranged parallel to each other, as shown. in appendix 5. 
Since each factor score in each component was associated. with 
a particular sediment sample, the location of each sample on 


each component was. found. Each component was then divided 


ne | oe Se 
4 y ‘SF 1 hth Me ete 

; ue : ' 

a | Fg - 

Pats Berobro siness sts tpo0x i> .~Yiledes _ nein igual Mod pai mie 

a 7 ie : ; ' a 

 Selistsb oxom ¢ 7 } #epeSagA is} exatoew onto 

; ‘ ’ 

Ca x I ok , 

a ainys otacwiah laefoatsg to: okt 

a; i 


j a } on : ena >a 4 or 
7 ark oa : a b' S x 
> LP w= 7 
; fig 2] satan eteevat ore ALeyi 
| : sam nokwelerxsn edt « 
we 7 rv pee oy 
. 
elds a to eSxpeb sis bas - aeidek 
é ‘3 a a 
Pe, > tacit ons. 93 i 
S56 Pits cau leva pio arte» f 
il ; i bs / ; ay ' a vi of 
Sr. tor- ) > ; Srienogmeos ipsa yd bee whe 
7 pe 
Bris * Ln ‘ Mt. OES 22 1S AOGMOS A ae 
Gsi2 BL) site a7] qT ,620uesert ey. t3nde39 xd 3%. 
Ae (picA pea Sastp-9 neuipmeienih sino 20th i y 
7 % 
doso vd betudixgros epesitsav heats oe x8 30.Le we ed it & 
snekoamon A .tasoxsg Ol oF isube se Mens 193983p a tee: 
Zod sfc {a - } pc re 5 Bi ge Bit Dic Leto seeds, veniya = cz, 


; : os 
bis .equorip 102 prismol sav fovyair dase: Bx idz _— 73 
NS - ar" Vy, oi 

agit.) 2 een iramapnieiannee 


ateil sansdognios 


| orig ai E5809 


it INO09 ogni ‘ers pttonde to Boristsm onr of >: 


a19W “snanegmes . 


os rx ay biel ; 


, = 
. 


Paw 


e*, 
7 in Ph 
40) elanas na Bo 


‘a 


teq svi aa erepioagnas <4 


-saubevoag 
na\ebatheal! 103062 
88 sso. foas os i4 2 
ee tose ak pa 2e: 
3 saacedl eo elem 2 sae nomibe oa % 


‘ aesoqiug: ovisaaors gts 


R 
202 aidasqeo 
ve 


Bia 
pat wold 8 ‘4 


0» anny foo ent 
data 
’ 7 7 4 
we 
: 2s 


ae 


? 


33- 
into positive and negative portions according to the values 
of the factor loadings in. each list. Subsequently, the 
samples located in the positive portion of the 2nd to nth 
component lists (where n is the number of components accepted), 
were traced over the whole range of the first component. This 
was done on the assumption that the first component, which has 
the greatest explained variance and the highest eigenvalue, 
will produce the most interpretative results for its whole 
range of factor loadings. In the remaining components, the 
positive range of the factor loadings was used more as a means 
of adding information to that found. in the first component. 
Hence, to provide a more reliable range of factor scores in 
the 2nd to nth components, only the positive values were 
considered. This first phase of locating samples on the 
component lists will give an indication of the behaviour of 
the components in terms of the varying magnitude of the factor 
loading attributed to each sample in each component. 

The second phase of checking involved. identifying 
every sample in the first component in terms of the original 
diagnosis of sediment type, made in.the field for each. sample. 
Those samples in the first component which were found in the 
positive portion of the 2nd to nth component, were also 
identified in these other components, and in the same way. 

In the course of locating and marking each sample, any 
distinctive groups or combinations of samples found within 
each component, and any groups or combinations identifiable 


in more than one component, were noted. The final grouping 
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of samples is then plotted in columns, component by. component 
not by marking individual samples but ina simplified form 
where groups of samples of the same sediment.type. are shown - 

as one unit in the column representing the component in which 
they were found (as shown in figures 10, 11 and.12). The 
clarity of the groups will give an indication of the degree. 

of separation achieved by this form of analysis: 

The final step is to isolate those samples which do 
not appear to fall into the same pattern as others with the 
same original sediment type diagnosis. In these cases field 
notes, and the original grain size data are referred to in 
order to identify a reason for the apparent anomalous 


behaviour. 


Results 

The initial inspection of the plots of the percentage, 
cumulative percentage, and logarithmic cumulative percentage 
grain size data, supported an earlier assumption that the 
grain size curves would be.too open-ended to permit use of 
Folk and Ward graphic statistics or moment measures. (See 
figure 8 for example). Therefore, the results of- the analysis 


outlined in the previous section (under the heading of method) 


are presented. 


Rank Order Listing 
The positions of the five sediment types, marked by 
individual. sample cases on the rank order list, are shown by 
sediment groups in figure 9. The length of the columns in 


figure 9 is proportional to the total number of samples 
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This shows the 75 80 and 90 
percentile ranges of occurence of 
the sample cases of each sediment 


type in the rank order list. 
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analysed. Each sediment type has been isolated from the 
complete rank order list so that the range of percentages 
over which each type is spread can be more easily seen, [In 
all cases there is bias of the median towards one end of the 
scale. The median of each type characteristically falls out- 
Side the 75 percent range of each of the other sediment types, 
except in the case of the fluvially derived sediments. The 
stability of these distributions is considered not good, 
especially where the number of observations is below forty. 
Interpretation of this grouping requires the generation of 
new groups. The new groups are based on suppositions. derived 
in the field, which have hitherto been. unmentioned. The first 
postulate concerns the overlap (b) of. the Laurentide till 
(Fig. 9) with the lacustrine types, and the overlap (a) of the 
lacustrine with the Laurentide till types (Fig. 9). The sup- 
position employed, is that mixing of the two sediment types 
has occurred. The samples found in the overlap in each of 
these sediment types are all from the same area of the basin, 
on the southwest flank of Tower Hill (Fig. 2), where possible 
mudflow conditions were identified in the field. At this 
point these two sediment types are in contact at the foot of 
a long steep slope. Furthermore, these samples are from two 
localities above and below the presumed contact of the lake 
sediments with the till. Hence, the overlap (b) of figure 9 
may be from the derivation zone of the mudflow on the. hill. 
slope above the lake sediments, whilst the lacustrine overlap 


(a) may represent the toe of the mudflow. The second 
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38 
supposition is that the samples from the Laurentide till 
overlap (c) with Cordilleran till samples which are together 
spatially, may in fact be Cordilleran till. 

In order to examine this new set of groups the Went- 
worth (1922) grain size classes of each sample were investi-- 
gated (in matrix MO-2). Observation showed that the medium 
to coarse sand and gravel classes were useful indicators of 
the distinctiveness of each sediment type in the Ree aude 
regions of the rank order listing. The lacustrine sediment 
samples in overlap (a) (Fig. 9) had a secondary low peaking 
in the medium to coarse sand and gravel classes with a magni- 
tude of 1.5 to 4 percent. This is similar to the Laurentide 
till samples, except that the peaking was much lower. for the 
samples in overlap (b) of this sediment. The Cordilleran 
tills show irregular peaking of the medium sands to gravels. 
and granules classes. The magnitude of percentages in the 
range medium sand to gravel. classes of the Cordilleran tills 
is slightly higher than in the same classes in most of. the 
Laurentide tills; except those samples in overlap (c) of the 
latter sediment type. It is therefore believed probable that 


the two suppositions referred to above are valid. 


Principle Component Analysis 


The correlation coeffecient matrix (Table 2) reflects 
the basic nature of the percentage values of the grain size 
data (matrix MO-1). As the percentage contribution of one 
grain size class tends toward 100 percent, so the other per- 


centages in the vector of observations will show a greater 
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40 
inverse relationship to this one class. The high percentage 
values of the first variable (clay-silt) indicates that the 
strong inverse relationships found between class‘ one and all 
Other classes (Table 2) are to be expected. The colinearity 
which may also pave weet oa in grain size data is also evident. 
in the correlation matrix, as adjacent values in Table 2 are 
mostly very similar, and the n to ntl classes are all highly 
correlated. The exception is found in classes 2 to 5 inclu- 
Sive. In this case zero values are common which suggests an 
independence of these classes to classes 6 to 15 inclusive. 

In classes 6 to 15 inclusive the original data values are. 
small and zeros are frequent, which would explain the high n 
by ntl correlation values of these classes once the first 
class has been considered. 

The eigenvalues. listed below (Table 3) support the 
acceptance of only the first four variables at the arbitrary 
level of 1.0 suggested by King, L. J. (1969). 

In terms of a percentage contribution of each variable 
or component to the total variance, an arbitrary level. of 
contribution of ten percent would reduce the acceptable number 
of components to three (Table 3). In fact this was the final 
number of components accepted. In the arguments presented 
previously using the rank order listing and the Wentworth 
groups, it was found that the third to sixth groups ef this 
scale were used to qualify fairly obvious distinctions rather: 
than absolutely determine them. Hence, the first and second 


Wentworth groups of clay-silt and fine sand respectively, 
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TABLE 3 


EXPLAINED VARIANCE AND EIGENVALUES 


FOR THE ANALYSIS OF ALL CASES 


Percentage Variance 


Component Eigenvalue Explained 
ip 6.327 42 
2 3.716 25 
3 25181 14 
4 1.233 9 
5 0.644 4 
6 0.360 2 
7 0.270 é 
8 0.130 1 
9 0.053 1 

10 0.046 2 
11 0.018 2 
12 0.013 ° 
13 0.009 
14 0.002 p 
15 0.000 v 
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will reflect the behaviour of all groups. These two groups, 
comprised of the first four grain size classes, or variables, 
are accounted for in the first three components of the prin- 
Cipal component factor listings. 

The factor loadings were arrayed in rank order with 
each component listing adjacent one to another (Fig. 10). 
When the first three components had been accepted, an arbi- 
trary datum or level of significance was established, using 
only values greater than unity. The first step of the inter- 
pretation involved checking the positive values of component. 
one with those of components two, three and four. Component 
four was included to-ensure that the major clusters were all 
accounted for. Those cases with positive values in all three 
components were found to have low percentages of clay-silt. 
In components two and three, all the fluvial samples had 
factor loadings greater than 1.5. This meant. that other 
sediment types could not form distinct groups in the greater 
than 1.0 range of the factor loadings. Thus, this first step 
was. repeated for the whole range of component one, against 
the positive ranges of components two and three. The desig- 
nated sediment type for each sample was also listed against 
the respective factor loadings in each component. The reason 
for this repetition is based on the tendency for the develop- 
ment of inverse relationships between class one (clay~-silt) 
and all others (see Table 2). In addition, the groupings of 
higher values in variable one will tend to produce lower 


ranges of factor loadings, an observation made by Gould (1967). 
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Therefore, investigation of groupings in the lower: range of 
component one may be justified. 

The results of this second phase of checking are now 
summarized cf figure 10. There are two dominant clusters. 
The first cluster is that of the fluvial, and ice wedge 
samples,.in the greater than unity range for ail components. 
The second cluster includes all the lake samples in the range 
less than zero of component one, and the range.1.1 to zero of 
component three. In component two only those lacustrine 
samples with values -1.6 to zero in component.one are found. 
This would include less than 25 percent of all the lacustrine 
samples and specifically includes those from the overlap (a) 
(Fig. 9) in the rank order listing with Laurentide till 
samples (from sites located in Area B, Fig. 2). In component. 
one these same Laurentide till samples (from Area B, Fig.2) 
range from -1.2 to +1.0, whilst in component two these till. 
samples overlap with lacustrine samples again. The Cordilleran. 
tills are found in the range +0.35 to +12.73 G6£ component one, 
except five which overlap with Laurentide samples in the range 
+0.35 to +1.0. This overlap is also found in. component three 
involving the same samples in each case. No Cordilleran till 
samples are positively ranked in the second component. Hence, 
these overlapping till samples (in.the component lists) lend 
support to the second supposition made earlier - that the 


Laurentide till of this overlap, which includes the same 
samples in the overlaps (c) and (d) of the two tills in the 
rank order listing (from sites in Area A, Fig. 2), may in 


fact be Cordilleran in type. Similarly, the overlap of 
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lake sediments from components one and two supports the first. 
supposition of mudflow action between Laurentide till and 
lake sediments in the vicinity of Area B (Fig. 2), 

The above procedure, involving principal component . 
analysis, was repeated for the two major overlapping groups. 
The first overlap is between the two till types and the 
second is between lacustrine and Laurentide till. 

The result of testing the Cordilleran. and Laurentide 
tills together yielded very distinct cluster centered on the 
zero value of the factor loadings of component one (Fig. 11). 
The Cordilleran tills have a characteristic positive loading. 
All the Laurentide tills have negative values in the factor 
loadings except for those taken from sites in Area A (Fig. 2) 
near Hoodoo Ridge and Tower Ridge. Those Cordilleran till 
samples with negative values were collected from north of 
Bragg Creek. Since these two localities have no samples 
taken between them there is doubt concerning the correct 
identity of the till samples from these two areas. These’ 
results are based solely on component one, which has-a 59 
percent explained total variance. Despite the fact that 
components two and three should be considered on the basis 
of the explained variance and associated eigenvalues, 
interpretation of each of these components reflected increased 
mixing of types within the factor loading lists. Since 


component one had a high explained variance and eigenvalue, 


(Table 4) interpretation was based on this alone. The 
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Component | Component 2 Fig ure ] ] 


GROUPED TILL SEDIMENTS 
ON LISTING FROM 


PRINCIPAL COMPONENT 
ANALYSIS 


Numbers indicate factor loading values at that point 


Gree! Cordilleran Till 


oo Laurentide Till 


Negative values 


Number of cases in each group. 


Cordilleran Till 37 


Laurentide Till 43 


Length of column is proportional to 


number of cases in each group. 
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TABLE 4 


EXPLAINED VARIANCE AND EIGENVALUES 


FOR THE TWO TILLS ANALYSIS 


Percentage Variance 


Component Eigenvalues Explained 
1 8.90 39 
2 3.43 23 
3 Pes 8 
4 0.69 4 
5 0.38 3 
6 0.20 1 
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correlation matrix for this phase of the analysis showed that- 
the first variable (the clay-silt fraction) ,- was. highly- 
inversely related to all other variables with values greater 
than 0.5. All other variables were related directly with the 
n to ntl and n+2 cases, all of which had coefficients greater 
than, 0.5. 

The results for the lacustrine and Laurentide till 
samples together (Table 5) yielded no such clear clusters. 
However, the results did produce certain useful characteris- 
tics. The associated correlation matrix had lower values 
than above. The relationships were similar overall, with 
locally high correlations in the low percentage variables. 
The eigenvalues and explained variances in each variable 
showed that the first four variables could be considered 
arbitrarily to be significant. The explained variances are 
very similar to those of the original analysis using principal : 
components (Table 3). This would suggest that these two sedi- 
ment types when combined, generated most of the variance of 
the. whole range of sediments used originally. These two sedi- 
ment types, which account for 70 percent of all the cases 
analysed, seem to create a bias in the variances. generated 
for each component. This bias is emphasized by the previous 
analysis, where when the number of cases of each sediment. 
type was similar, the ranges of values of the first component. 
were different (Figs. 11 and 12), and distinct clusters were _ 
generated. In the present case the first and third components 


show similar and reasonably strong grouping (Fig. 12). The: 
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TABLE 5 


EXPLAINED VARIANCE AND EIGENVALUES FOR THE: 


LACUSTRINE-LAURENTIDE TILL ANALYSIS 


Percentage Variance 


Component Eigenvalue Explained 
i 6.45 43 
2 3.45 23 
3 1.88 13 
4 ig Ei 8 
5 04 72 ° 
6 0.50 3 
7 0.35 3 
g 0.20 1 
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Indicates the same case 
located in each component 


listing 


Figure 12 


GROUPED LACUSTRINE AND LAURENTIDE TILL ON 
LISTING FROM PRINCIPAL COMPONENT ANALYSIS 
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second and fourth components show similar but mixed groups. 
In both cases the interpretation of the first principal compo- 
nent analysis were not improved upon. In general therefore, 
this group probably reflects the bias produced by the number 
of cases. In addition, within group variances are. important, 
but not great enough to produce distinct clusters within a 
high initial total variance. Hence, the high percentage 
values of the first variable (the clay-silt fraction) of the 
original data is insufficient to produce good differentiation 
of these two sediment types. 

The final allocation of. sediment type was made by 
continued investigation of the overlapping cases; and field 
notes and grain size data were also used. The final results 


were plotted on the surficial geology map (Fig. 21). 


Conclusions 

This experiment has shown that the behaviour, of the 
three components is as follows. The first. component tends to 
produce a range of factor loadings whose values are inverse 
to the class one grain size percentages. Only percentages of 
the clay-silt class less than 40 percent produced positive 
factor loadings. The second component is more sensitive to. 
larger variances in the second and possibly third classes of 
grain sizes. This continues to produce inverse factor 
loadings but emphasizes stronger variance in class two with 
respect to class one, and thereby eliminates Cordilleran tills 
and most lake sediments. The third component, which brings 


fluvial and lacustrine samples into juxtaposition, appears to 
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group those cases which have extreme values in the clay-silt. 
percentage, thereby emphasizing the maximum variance of vari- 
able one. The fourth component is far more varied in behav- 
iour and is far less easy to interpret. This component has 
put ali the tilis and lake samples into small clusters with 
no sorting between types. Hence, the analysis was discon- 
tinued in the fourth component. 

The analysis, which was. used primarily for interpre- 
tative purposes, has shown that discrimination based on prin- 
cipal component analysis between a number of sediment types 
is possible. In addition, it indicates those cases where 
there is doubt as to the correct sediment type allocation and 
provides sufficient information for the reasons causing this 
doubt to be identified. In this study the analysis indicated 
that eleven samples had been wrongly identified, and that 23 
others required additional information before definite identi- 


fication could be made. 


Summary 

A model for the testing of hypotheses relating to the 
determination of sediment type and the spatial significance 
and pattern of the sediments could be built upon the approach 
used in this study. Many workers have used factor loadings 
as data for the statisitical mapping of phenomena. A similar 
approach was used by Kendall (1939) for crop productivity in 
England. In this application it is felt that the shape of the 
area studied and the density and distribution of points, do 


not provide sufficient. grounds for applying techniques. of 
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spatial analysis such as canonical trend surface analysis 
(Lee, 1969) to the data obtained. Sampling methods which 
would however, permit such research are as follows: 

1. A regular or more even distribution of sample 
points for any given scale. The juxtaposed nature of con-. 
trasting sediments encountered in-this study suggest that a 
small area be used initially, so that mapping.criteria (e.g., 
for isoline) may be. ewaluated. (Ogawara and Yamazaki, 1952). 

2. More complete granulometric data will give an 
improved grouping of factor loadings within components and 
allow the mapping of component by component. 

3. The inclusion of additional: variables. 

4, Finally, the researcher must develop criteria by 
which the delimitation of areas may be achieved according to 
basic geomorphological considerations. 

A possible avenue for future research could involve 
two or three new variables. The object would be to identify 
spatial patterns similar to those found by conventional 
methods. 

In conclusion therefore, the finding is that principal 
component analysis can be used to discriminate between known 
sediment types. It may also be used to identify sediment 
groups not easily distinguishable. The usefulness of-prin- 
cipal component analysis lies in its flexibility in terms of 
identifying factor loadings for every sample. The differences 
in grouping within each component are the expression of within 
group variances and hence this technique conveniently expresses 
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the limitations of the original data. The major disadvantages 
of principal component analysis seem to be the difficulty of 
applying significance tests, and the more fundamental fact 
that principal component analysis does not give any idea of 


error variance except by interpretation. 
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CHAPTER III 


RECENT GEOLOGICAL HISTORY OF THE ELBOW 


RIVER BASIN. 


Several general statements have been made by Beach 

(1943) concerning the study area, and the author concurs with 
these. The first emphasises structural control. 

"the markedly different resistance to erosion offered by 

the various lithological units of the bedrock has been 

the controlling factor determining the nature of the 

topography .. . The ridges mark the outcrops of massive 

sandstone and conglomerate beds and the flat-bottomed 


valleys those. of formations. composed predominately of 
shales." 


Beach also recognised a major discordant feature. 
"The foothills part of the region is dominated by one 
major topographic feature, Moose Mountain Dome. .. . 
Streams have cut several deep canyons across the dome." 

Except for a few Tertiary and late Cretaceous outcrops. 
found east of Last Break Hill, these quotations, summarise the 
major morphological units of the area. The Tertiary strata 
which underlie this area of lowest relief (termed Prairie in 
Chapter I) are soft sandstones and shales of the Edmonton, 
Upper Alberta, Belly River and Paskapoo formations. 

Figure 13 illustrates the dominant northwest to. 
southeast trend of the geological strata and structural 
geology. The highly folded and faulted foothills regions 
bring the Palaeozoic sandstones, dolomites, dolomitic lime- 


stones and conglomerates into juxaposition with Mesozoic 


shales and sandstones. of a less resistant nature. To the. 
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east, the late Cretaceous rocks are faulted where they overlie 
the eastward extending Mesozoic and Palaeozoic rocks. The 


Elbow River is discordant with the geological trends 
The Problem 


This discordance between the course of the Elbow River 
and the orientation of the geological structure within the 
foothills poses the problem of how this situation arose. The 
four solutions considered are outlined below. 

Ll. The concept of river capture was considered by 
Sissons (1960) to be important in the evolution of drainage 
in southern Scotland. In the study area, river capture may 
be responsible for the discordance of drainage and structure. 
This process could facilitate diversion of flow from struc- 
turally controlled valleys into valleys which are discordant 
to structure. In order to test the hypothesis of river 
capture, alignment of present rivers and the presence of 
Windgaps should be investigated. 

2. The breaching of watersheds by glacial ice is an 
alternative method by which rivers may be diverted from one 
drainage system or channel to another. Similar features to 


those associated with river capture together with glacial 


features will identify this type of process. 
3. A process frequently considered in a region in 
which there is a problem of discordance of river to structure, 


is superimposition. The application of the superimposition 
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concept to explain discordance is described by Wooldridge and 
Linton (1955). The concept implies the existence of a drain- 
age pattern on a surface prior to the initiation of rejuvena- 
tion of the drainage system which is superimposed. This 
surface may be either a re-cycled peneplain or a recently 
emerged sea-floor. The passage of time sees the down cutting 
of river channels due to isostatic and/or eustatic change, 
such that the rivers will eventually be. established in rocks 
of an older period than those from which it was initiated and 
discordant with the structure of these rocks. At this point 
the river is said to be superimposed. This process specifi- 
cally excludes orogenic effects but admits epeirogenesis. The 
evidence cited for the process may include remnants of the 
original surface or sedimentary cover from which the drainage 
was developed. Without information relating to the presence 
of the original surface, it becomes very difficult to infer 
that this process has operated within a region. 

4, The process of antecedence implies that a river 
is able to maintain its course whilst the bedrock. and the 
ground surface over which it flows are subject to orogenic. 
disturbances, The river is able to downcut at a rate which 
is greater than the rate at which structurally controlled 
ridges are formed. The process has been described by Sparks 


(1964) as " the last resort of a geomorphologist seeking 
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to explain inconsequent drainage, .« - ." because it involves 


a guesses as to the rate of uprising folds and the rate 


of downcutting of rivers, about neither of which is very much 
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known." However, if the period of formation of geologic 
structures can be associated with the period which post-dates 
the beginning of the present phase of sub-aerial denudation, 
then antecedence may be implied. In this situation the 
presence of bedrock benches along the valley sides would help 
to substantiate the hypothesis of antecedence. 

The background of each of these four possible solu- 
tions to the discordance of rivers is geological. The deci- 
Sion that any one particular process is responsible for the 
Situation will depend upon the morphological information 
available. The geological and morphological factors which 


have been found in the study area follow. 


Late Mesozoic and Tertiary Geology of the Study 
Area and Environs 


Late Cretaceous 

The basal formations of the late Cretaceous geology 
show evidence of a marine environment. Above the basal 
formations are sandstones and conglomerates in which Beach 
(1943) found pieces of bone (at Quirk Creek, in the Elbow 
basin) later identified as dinasaur bones (a femur and jaw 


bone). The sandstones form a part of the Edmonton formation. 


Tertiary 


The first well established event in the Tertiary era 
is the depositional phase of the Paskapoo formation. Accord- 
ing to Tozer (1956) the Paskapoo formation is of sub-aerial 


origin throughout in the vicinity of the Elbow River, because 
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the recession of the marine environment had begun in late 
Cretaceous times as the Wapiabi Sea became restricted to the 
Alberta syncline. Formation of the Paskapoo deposit continued 
through Paleocene times into the Eocene period. Thereafter 
deposition was irregular as the earth movements of the 
Laramide orogeny, which began in mid-Eocene times, interrupted 
the spatial continuity of the Paskapoo evident in Paleocene 
times. The Laramide orogeny continued throughout the rest of 


the Tertiary period. 


Distribution 

The map of the geology of the Elbow basin (Fig. 13) 
shows that the Paskapoo formation is found in the Prairie 
region alone. The Edmonton formation, which includes the 
Bearpaw equivalents, is found along the boundary zone between 
the foothills and the Prairies (Beach, 1943). The Belly River 
formation outcrops in a synclinal structure on the west side 
of the Edmonton formation (Fig 13). The Belly River formation 
is also found faulted against Palaezoic rocks of the Rocky 
Mountains, along the Rocky Mountain-foothills boundary (Beach, 
1943). All rocks older than the Belly River formation outcrop 
west of Bragg Creek. 

In the mapped area of figure 13 are seven major struc- 


tural units. To the west are the Palaeozoic strata of the 


Rocky Mountains (Z), (Fig. 13). Highly faulted and intensely 


folded, this area has high relief and a complex structure. 


To the east is an area (Y) separated from the Rocky Mountains 


by a major thrust, known as the Nihahi Thrust because of its 
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location on the eastern flanks of Nihahi Ridge. This is the 
area of the Forgetmenot Anticline (Fig. 13). North and east 
of the latter is the prominent structural and topographic 

form of the Moose Mountain Dome (X). Both of these anticlinal . 
structures cause outcropping of Palaeozoic strata (Fig. 13). 
East of Moose Mountain Dome is a faulted syncline in late 
Mesozoic strata (W). This area is flanked on the eastern 

Side by highly faulted strata of early-middle Mesozoic rocks, 
which form a low anticlinal structure (V). The eastern limb 
of the latter slopes eastward, and is covered by late Mesozoic 
strata which are highly faulted in area (U). The most easterly 
area (T) is not folded and, except for the juncture with area 
(U), not affected by faulting. The area (T) is formed by the 


Paskapoo formation. 


Geological History of the Environs of the Study Area 

The work of McCrossan and Glaister (1964) shows that 
during the period in which the Belly River formation was 
being deposited, the area of the associated marine environment 
was decreasing (Fig. 14, A-F). By the end of this period, 
the foothills of Alberta were exposed to a sub-aerial environ- 
ment, whilst the marine basin extended south from south- 
eastern Alberta. The subsequent depositional phase associated 
with the Edmonton formation was marked by a western shoreline 
in a similar position to that of the end of the Belly River 
period (Fig. 14), but the marine basin now extended much 
further east into southern Saskatchewan and southwestern 


Manitoba. The western margin of this basin, which was located 
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UPPER CRETACEOUS PALAEOGRAPHY 


OF WESTERN CANADA 


FIGURE 14A 
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FIGURE 14C 
MILK RIVER SEA ( Lower Campanian ) 
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LOWER BELLY RIVER SEA ( Middle Campanian ) 
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FIGURE 14E. 


UPPER BELLY RIVER SEA ( Middle Campanian ) 
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Source : McCrosson and Glaister , 1964, p. 187. 
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FIGURE 14F 
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in the study area, oscillated locally such that marine deposits: 
are found overlayen by sandstones and conglomerates, as Beach 
(1943) describes in the Elbow basin. Russell (1957) suggested 
that the oscillations were caused by epeirogenetic adjustments 
of the land surface. Hence, at the end of the late Cretaceous. 
period and the beginning of the Tertiary era uplift associated 
with the pre-Laramide disturbances began, and marks the exten- 
Sion of the sub-aerial environment already existing in eastern 
British Columbia and western Alberta over the larger part of 
southern Alberta and specifically, the whole of the study 
Besa 

The Paleocene period is suggested by Russell (1932, 
1950 and 1957), Elliot (1958), and Tozer (1956) as the period 
of formation of the Paskapoo stratigraphic unit. The Laramide 
orogeny began in Eocene times and continued throughout the 
Oligocene period. No distinct depositional phase is associated 
with these two periods. Minor earth movements are thought to 
have occurred in Miocene and Pliocene times by McCrossan and 
Glaister (1964). In addition there is evidence for sub-aerial 
depositional phases during these two periods for several very 
local Br can To the west of Calgary (30 miles) are the 
Crawford Hills which are capped by gravels (at an elevation 
of 4,400 feet) as are the Hand and Wintering Hills of central 
Alberta. The problems associated with the interpretation of 
gravels found on isolated upland surfaces in Alberta have been 
widely discussed by Westgate (1965) and Crickmay (1965) for 


the Cypress Hillis in southeastern. Alberta. 
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Topography of the Study Area and Environs 


The topography of the Elbow basin shows a general 
decline in the elevation of summit peaks from west to east 
(Pig, 15). There are distinct areas of homogeneous topog- 
raphy (Fig. 15) close to an elevation of 4500 feet, and in 
the eastern foothills several flat topped hills rising up to 
5000 and 5500 feet were found. No bench features were found 
associated with the bedrock topography. 

The bedrock topography and pre-glacial channels of 
the southern part of Alberta have been mapped by Geiger 
(1965), Farvolden (1963), Stalker (1961) and Carlson (1969, 
1969A and 1970). None of these works have mapped a pre- 
glacial valley in the present Elbow basin. However, Farvolden 
(1963) was- able to trace the pre-glacial valley of the High- 
wood River (in southern Alberta, 50 miles south of Bragg Creek) 
from the foothills onto the Prairies, although he points out 
that there is a break in the alignment of the preglacial valley 
between these two regions. In the case of the Elbow River, 
Farvolden traced a bedrock valley westward from a point seven 
miles south of Calgary up to the beginning of the foothills 
near Priddis, but not into the foothills. In the work of 
Geiger (1965), Farvolden (1963), and Stalker (1961) pre- 
glacial bedrock channels, other than those of the Oldman, Bow 
and Red Deer rivers, are not traceable into the foothills. 

There is a lack of additional morphological informa- 
tion for the Elbow basin in the literature. In view of this, 


analysis of the topography was carried out in order to 
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identify morphological relationships between high and low 
elevations, and the major topographic units. There are 
several techniques which may be used for this purpose, many 
of which have been discussed by Pannekoek (1967). In this 
study two techniques are used, namely isoline mapping and 
trend.surface analysis (King, C.A.M:, 1969). The data con- 
Sists of the elevations of summit peaks from the Elbow basin. 
The method of collection of this data is given in appendix 6, 
part. 13 

It is generally hypothesised that summit peaks are 
remnants of an older surface on which drainage was. initated 
(an erosion surface). Based on this assumption many workers 
in various parts of the world have used summit peak informa- 
tion in order to identify an erosion surface. In this study 
the peaks are assumed to be useful indicators of an older 
phase of drainage initiation. It is considered that the 
distribution of summit peaks and the relationship of one peak 
to another will be more reliable for interpretative purposes 
than the specific location and elevation of a single peak. 

The two forms of analysis to be used both. establish 
a form of continuity between data points. It should therefore 
be stated that the contour lines drawn are generalisations of 
an unknown morphological form, and that the interpretations 


based on the results are therefore made about the distribu- 


tional aspect of the summit peaks. 


Isoline Map of Summit Peaks 


The data is plotted (as in figure 16), and a system 
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OF REFERENCE POINTS FOR 


FIGURES 16 and 18 
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of contouring applied to generate a best fit for the data 
points obtained. The location of contours is made subjec- 
tively, according to the best fit generated by a given set 


of circumstances. 


Trend Surface Analysis of Summit Peaks 

An account of the statisitcal procedure for trend 
surface is given in appendix 6, part 2 (see King, C.A.M., 
1967). The computations were made using a computer program 
prepared by O'Leary, Lippert and Spitz (1966). 

The assumptions made when applying trend surface 
analysis to summit peaks in the Elbow basin are listed below. 

(ary) The summit peaks will reflect the basic spatial 
pattern of high and low elevations of the area, despite 
lowering of the points since the initial phase of incision. 

(ii) Since no bedrock benches were found in the 
basin, the surface fitted by trend surface analysis will not 
reflect any one geomorphic surface (Brown, 1957), but an 
average "trend" surface which may be used to interpret the 
major valley trends. Hence, the surface derived is not the 
erosion surface, but an indicator of the form of the surface 
being considered. 

(iii) The nature of the "trend" surface further 
emphasises that the trend is not the surface of initiation 
of erosion, which was conceived as a multitude of phases of 
development dependent on uplift and most probably reflected 


in benches 


(iv) The asset of the trend surface is. that 'noise' 
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70 
(Board, 1965) or local components (Allen and Krumbein, 1962) 
of the topography are removed. These are presented as 
residuals. 

(v) The residuals are assumed to reflect structures 
which will either indicate the main reasons underlying the 
discordance of the river to the geological structures, or the 
independence of topography from geological structure. 

The reasons which necessitate the above assumptions 
were most clearly described by Tarrant (1970) in his criticism 
of. the work of C.A.M. King (1969). 

In 1969 King applied the technique of trend surface 
analysis to summit peak data for the Alston and Askrigg blocks 
of the Central Pennines of England. Tarrant (1970) indicated 
that the surface fitted by trend surface analysis does not 
envelope the summit peaks, as King suggested it would. 
Instead, the surface is fitted in the same way as a least 
squares regression line. The statistics of trend surface 
analysis require that the surface be arranged such that the 
magnitude of the deviations below the surface equal those 
occurring above the surface. The positive residuals are not 
identifiable as monadnocks. Therefore, the trend surface does 
not reflect the form of the erosion surface as C.A.M. King 
(1969) had implied. However, in trend surface analysis the 
relationship iS such that more than two observations will 
affect the relationship of a third and locationally inter- 
mediary point. In addition, the generalising effect of trend 


surface analysis suggests that forms evident within the trend 
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of the surface are dependent upon the data and not the tech- 
nique. This does not mean that details of the elevation of 
the trend surface specifically identify the location of the 
components of the erosion surface, nor does it indicate that 
an erosion surface existed. It does however, identify any 


major topographic phenomena in the area studied. 


Results of Analysis of Summit Peaks 


The Isoline Map 

The results are shown in figure 16. The major 
features are indicated by letters (the names of which are 
iMuoxed on the page following figure 16) 

The pevise shown on the isoline map (Fig. 16) 
increases westwards from a line north-south through Bragg 
Creek. To the east of this line the relief is low and reason- 
ably smooth. This may indicate that statistically the data 
for peaks east of this line reflect colinearity to a greater 
extent than data for peaks west of the line. 

The topographic forms found west of the north-south 
line through Bragg Creek are largely related to geological 
structures, for instance (Fig. 16) Moose Mountain Anticline 
(G), Forgetmenot Anticline (F), Nihahi Syncline (R), Mount 
Howard complex (Fairholme Anticline (E)), Kananaskis fault 
line (C), and the Barwell Syncline (S). The topography of 
the Rocky Mountain-foothills area would seem therefore, to 


be dependent upon geological structure. 


The topography east of the north-south Bragg Creek 
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72 - 
line cannot be closely related to geological structure, and 
may be considered to be independent of structure. 

Preliminary investigation suggests that the Elbow 
River is closely aligned to topographic depressions, or 
valleys in figure 16. Further inspection shows that at only 
One point does the river conform with such a depression. 
This point is marked (X) on figure 16. Elsewhere, the river 
is aligned to, but not centred on the low areas proximal to 
its course. At point (Y), the Elbow River. appears to diverge 
from a major valley (P) which can be traced eastwards to 
Priddis (Fig 16). The Elbow River flows northeast from point 
(Y) and cuts across the entrance to Bragg Creek valley (Q). 

The rivers east of the north-south line through Bragg 
Creek are accordant with the valley forms mapped from the 
summit peaks. It is noticeable that the present Priddis 
Creek occupies the major valley (P), but not completely. As 
with other present day rivers in the eastern area, Priddis 


Creek headwaters are all in the same area. 


The Trend Surfaces 

The results of the analysis are shown in figures 17 
and 18, where only the sixth order surface and associated 
residuals have been mapped. Inspection of the first through 
sixth order surfaces showed that the sixth order surface was 
the most useful for the reasons given below. 

(1) The topography as shown in figure 15, suggests 
that spatial autocorrelation is present. This means that 


simple trends of the summit peaks will yield no more 
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information about highs and lows of the data than that which 
is already evident, namely increasingly higher values west- 
ward from Calgary. Inspection of the 1 through 6 order. sur- 
faces showed that only the fourth, fifth and sixth show any 
recognisable patterns of high and low areas within the general 
east-west trend. 

(2) The high correlation coefficients for all surfaces 
(Table 6) further supports the assumption of autocorrelation. 
However, the variation explained, the standard deviation, and 
the coefficient of determination for the fourth, fifth and 
sixth order surfaces are all quite close. These statistics 
for the fourth and fifth order surfaces are very close, whilst 
the sixth order statistics show a recognisable difference. 

(3) All of the surfaces were found to be significant 
at the greater than 99.9% confidence level. This is thought 
to be attributed to the large numbers generated for the total 
variation. The statistics used to find the significance level 
is the 'F' test, (see appendix 6) where.the larger mean square 
value is divided by the smaller mean square value. In this 
case if the mean square values are generated from large varia- 
tion values and there is a large difference between the two 
sets of degrees of freedom, then since for each surface the 
numbers of degrees of freedom are fixed, the magnitude of the 
number of observations will determine how sensitive this form 
of significance test will be for the analysis in question. 


In this study the F value for the sixth order surface was 


found to be 544.0, and the tabled F value was 1.74 for 30 and 
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400 degrees of freedom, at the 99.9% level. 

(4) The first order surface is accepted as being 
indicative of a general, regional component in trend surface 
analysis. In this study the first through fourth order sur- 
faces show similar trends, but very little variation within 
the surfaces. The fifth order surface contains troughs, 
whilst stiil maintaining the trend of the surface. ‘In the 
Sixth order surface these ivefighkearc more distinctive and 
the general trend of the data. is maintained. Since trend and 
patterns of high and lows of the peaks within the surface are 
the factors originally assumed to be useful in this study, 
the sixth order was used to interpret forms evident within 


the surface. 


The Sixth Order Trend Surface 

The sixth order surface (Fig. 17) shows. increasing 
elevations westward. The strike of the surface in the 
central foothills portion of the map is approximately north- 
west to southeast, and the dip is oriented southwest to 
northeast. The Elbow River is approximately aligned with the 
dip, though local changes in direction away from the dip 
occur. In the eastern part of the map (Fig. 17) the surface 
indicates the presence of a major valley in this area which 
begins close to Priddis and extends eastwards. There is a 
feature which resembles a col in the extreme west of the map 
(Fig. 17). This col extends westward from Mount Howard, and 
is therefore located well within the Rocky Mountains. As the 


trend surface represents major regional trends, according to 
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Allen and Krumbein (1962), it is assumed that the feature shown 
in figure 17 are regional rather than local components of the 
topography. Hence, the central portion of the map is indica- 
tive of the regional topographic components of the foothills. 
In the case of the Rocky Mountains and the Prairies the size 
of the sample area is too small in comparison to the dimen- 
Sions of these two regions, and the trend surface may not 


therefore be representative Of these areas. 


The Sixth Order Residuals 

The pattern of the residuals plotted and isolined in 
figure 18, are similar to the isoline map of the raw data 
(Fig. 16). The letters used to locate features in figure 18, 
are Similarly shown on figure 16. Comparison of these loca- 
tions will reflect the degree of similarity of the two maps. 

To the west of the line north-south through Bragg 
Creek, the following features are shown similarily on each 
map (Fig. 16 and 18). The high areas A,B,D, E,F and G. The 
low areas C,R,S and P. The residual plot (Fig. 18) shows a 
Similar pattern of highs and lows west of a line north-south 
through Moose Mountain (G), to that of the isoline map (Fig. 
16). To the east of this line the pattern is different 
between the two maps. The relief is lower east of this line 
on both maps. On the residual plot there is shown a long low 
ridge trending NNW. to SSE. through Bragg Creek (Big. 228) 
which is not found in the isoline map (Fig. 16). To the east 
of this ridge is the high area marked (K). Between these two 


ridges there are no data points, which means that there 
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were no clearly identifiable peaks in that area. Hence, the 
isoline map which tends to produce continuity between points 
of equal value, shows the area between the ridges as an homo- 
geneous flat area, whilst the effect of the residual plot was 
to bridge the gap as in a regression analysis. However, both 
maps may represent this area equally correctly, as both show 
a low area bounded by a ridge. 

To the east of the line north-south through Bragg 
Creek, the features on the residual map (Fig 18) are less 
clearly identifiable with the features in the same seen on 
the isoline map (Fig. 16). However, the basic forms are 
recognisable. 

The residual map tends to show that the Elbow River 
is more discordant with the topography of the summit peaks 
distribution than does the isoline map. The difference in 
behaviour of the two types of mapping used above, suggests. 
that this discordance on the residual map is exaggerated. 
Alternatively, the continuity of form produced by the residual 
map may indicate an underestimate of the degree of discordance 
shown on the isoline map. The difference between the discor- 
dance shown on each map is that on the residual plot the lines 
of the discordant features are at right. angles to the general 
line of the river. On the isoline map the discordance is one 
of the displacement between the valley and the river. Since 
for Priddis Creek the river and valley are accordant on both 


maps, two types of discordance of the Elbow River may exist. 
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Summary of the Topographic Analysis 

In terms of the interpretative qualities of each of 
the two techniques used, isoline mapping of the summit peaks 
was more reliable. The doubts concerning the trend surface 
analysis centres upon statistical fitting of the regressed 
surface and the degree of generalisation achieved. As relief 
decreases the behaviour of the polynomial fit appears to vary 
more than in the case of the isoline map. The trend surface 
analysis produces greater continuity of form in high relief 
areas, whilst the isoline map reflects smaller and more local 
features in Similar areas. 

There is reason to believe that an older surface of 
erosion may have existed, for the following reasons. 

The trend surface (Fig. 17) shows an accordance of 
the present Elbow River with the dip of a surface which has 
a. high degree of correlation with the summit peak data. 

The topography of the foothills is strongly influenced 
by geological structure. The Elbow River is discordant with 
geological structure, and also with the pre-Pleistocene 
topography in the foothills. To the west of the Elbow River, 
the author has noticed a greater accordance of rivers with 
structure, as for example in the Kananaskis valley. To the 
east of the foothills in the Elbow basin the topography is 
less dependent upon geological structure and the rivers show 
a greater accordance with the topography of the summit peaks. 

A point of explanation concerning the discordance of 


the Elbow River in the foothills downstream of point (Y) in 
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figure 16, is that below this point the Elbow valley has a 
thick surficial infill of gravels and till. The asymmetry 
of the river in the valley shown in figure 16 may. therefore 
be due to glacial interference. This point will be discussed 
further in chapter 4. Hence, the residual map may reflect 
the major form of discordance between the present river and 
the topography. 

A final point is that approximately fifteen miles 
west of the headwaters of the present Elbow River is the 


continental divide. 


Conclusions for the Topographic Analysis 


There was no sedimentary cover which extended uni- 
formly over the whole study area and its environs, and 
three basic structural and topographic units are evident. 


These are; 


The Rocky Mountains, east of the continental divide 
and west of Nihahi Ridge (Fig. 19, east of region l). 

The foothills, east of Nihahi Ridge and west of a 
line north-south through Bragg Creek (Fig. 19, regions 1 and 
2). 

The Alberta Prairie region, east of Bragg Creek (Fig. 
19, region 3). 


Interpretation of the Tertiary Geomorphic 
History of the Elbow Basin and Environs 


The end of the late Cretaceous period was marked by 


a marine recession from the vicinity of the study area. In 
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other local areas of Alberta marine sedimentation continued. 
Tozer (1956) described marine sediments at several localities 
along the Prairie-foothills margin, but noted that the contact 
between the Edmonton and Paskapoo formations showed a marked 
break in the Elbow basin, ". . . 10 miles south of: Cochrane 
> co e e The Paskapoo at. this basal level consists-of-". . 

- conglomeratic beds. .. . On Elbow River. a large bone frag- 
ment, presumably dinasaur, was collected 80 feet stratigraphi- 
cally above this conglomerate .. .." No shales have been 
found, which suggests that recession and subsequent sub-aerial 
deposition occurred in the Paleocene. Prior to this recession 
the late Cretaceous marine environment had not been extensive, 
and the western margin of the sea was in the vicinity of the 
foothills of Alberta (McCrossan and Glaister, 1964). Hence, 

a sub-aerial environment had existed in the area now known as 
the Rocky Mountains and the western edge of the Alberta foot- 
hills, prior to the Tertiary period... Thus, in Paleocene. 
times sub-aerial denudation and deposition extended over the 
study area. 

The Laramide orogeny began in mid-Eocene times 
(McCrossan and Glaister, 1964) and continued throughout the 
Oligocene period. Minor earth movements occurred in Miocene 
times and the orogeny culminated in one last phase of move- 
ment in Pliocene times. 

A sub-aerial environment persisted throughout the 
Eocene, Oligocene, Miocene and Pliocene periods. Westgate 


(1965), working in the Cypress Hills of southeast Alberta, 
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associated the upland surface of this area with a long period 
of sub-aerial denudation which lasted from mid-Paleocene 
times until the end of the Tertiary era. Crickmay (1965) 
postulated panplanation, a process of lateral planation by 
rivers, (Crickmay, 1933) as the cause of the flat surface. 
It was also suggested that by this process and late Tertiary 
rejuvenation the Cypress Hills became an isolated upland area. 
The Cypress Hills are on the eastern flank of the late 
Cretaceous~Paleocene "Alberta Syncline." The Elbow basin is 
on the western margin of this syncline. Therefore, it is 
expected that drainage flowed toward this synclinal area in 


southern-central Alberta. 


Conclusion: The Pre-Pleistocene Evolution 
of the Elbow River 

In the study area two pre-glacial valleys are identi- 
fied. One is in the Prairie region (Priddis Creek valley) 
and the other in the foothills (Elbow River valley). Neither 
of these valleys can be traced across the foothills/Prairie 
boundary. The presence of a pre-glacial valley and a 
Paleocene river in the foothills, both of which approximate 
the direction of the present valley of the Elbow River, sug~ 
gests that the Elbow River was antecedent to geological struc- 
tures formed during Eocene and Oligocene times. The present 
dependence of topography on geologic structure inothe: foot 
hills of the study area may further support this hypothesis. 
However, without more conclusive evidence within the study 


area of benches which may only be attributed to antecedence, 
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this hypothesis remains conjectural. 

An alternative explanation may be that of post- 
Laramide superimposition from a peneplain. In this case the 
assumption is that during Eocene and Oligocene times, 
topographic forms created by orogenic disturbances and pre- 
Laramide topography were reduced to a peneplain by the begin- 
ning of Miocene time. This peneplain may be reflected in the 
highly correlated summit peaks found in the trend surface 
(Fig. 17). Subsequent rejuvenation of drainage during Miocene 
and Pliocene times may have been caused by earth movements 
along the 'Disturbed Belt,' as mentioned by Tozer (1956). 
Hence, superimposition of late Oligocene/early Miocene drainage 
patterns may have taken place. A similar process is described 
by Westgate (1965) for the Cypress Hills, in SW. Alberta. 

In order to substantiate either one of these hypotheses 
two types of information are required. The first concerns 
more detailed investigation of surface and bedrock topographies 
in order to accept or reject the hypothesis of antecedence. If the 
hypothesis of antecedence is rejected, then a second source of 
information is required to substantiate the hypothesis of 
superimposition. This would involve a search for clastic 
sediments which may be attributed to the peneplanation stage 
of the Laramide geologic structures prior to possible super- 


imposition in Miocene and Pliocene times.. 
Post Tertiary Geomorphic History of the Elbow Basin 


The uncertainty of events within the Tertiary era of 
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Alberta appears to carry through into the Pliestocene era. 
Despite extensive glaciers over North America during the four 
major glaciations attributable to North America only the last, 
the Wisconsin, is discussed for Alberta with any certainty 
(Stalker, 1960). 

The pre-Wisconsin landscape was modified by Cordilleran 
and Laurentide ice masses and associated deposits. The extent 
of the modification and/or change which took place is not 
clear, since knowledge of the pre-Wisconsin landscape is vague. 
The amount of change due to Wisconsin ice alone is difficult 
to assess in the light of a possible earlier glaciation of the 


study area. 


Morphology and the Surficial Sediments 

The study is considered in three regions in this 

section. 
(i) West of McLean Hill. 
(ii) East of McLean Hill and within the foothill region. 
(iii) East of Elbow Bridge. 

(i) The first area (1B) is west of a north-south 
line through McLean Hill (Fig. 19). This area has the greater 
relief, and evidence of glacial erosion. In figure 20 the 
major morphological forms have been mapped. using field observa- 
tions and aerial photographs. In order to illustrate differ- 
ences in the morphology of the study area, the mapped area of 
figure 20 has been extended westward as far as Mount. Glasgow. 

A prominent morphological feature is the concentration 


of cirque forms in the higher relief areas close to and west 
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Of Moose Mountain (Region 1A and Westwards; Fig. 19). A 
second prominent feature is the extensive and high relief 
ridges which controlled ice movement. To the west of Elbow 
Falls the valleys are typically 'U' shaped with valley side 
moraines extending well above the valley floors. Within this 
area are three main cirque locations. (1) Forgetmenot moun- 
tain, (2) Mount Glasgow and west of Nihahi Ridge, and (3) 
around Moose and Jumpingpound Mountains (Fig. 20). Since ice 
probably did not spill over into the Elbow system from the 
Kananaskis valley (Pers. Comm. Mr. M. Walker, University of 
Calgary), the above were the source areas of Cordilleran ice 
which occupied the Elbow valley in Wisconsin times. The main 
directions of ice movement proposed are shown in figure 20. 
The movement of ice across the broad valley immediately east 
of Nihahi and Compression Ridge may be attributable to ice 
blocking at the north end by Bow and Kananaskis ice, as well 
as the pre-existence of river valleys which traversed east 
via Canyon and Prairie Creeks, and via the Elbow and Little 
Elbow valleys (Fig. 20). Therefore, the ice would probably 
have been thickest close to Elbow Falls, because at this 
point most of the ice from all sources was confined in a 
narrow valley. 

Bast of this’ point. (in region) 18%) Fig. 19) the main 
valley is broad and occupied by fluvial sediments varying 
from 50 to 150 feet thick. The valley walls consist of 
truncated spurs which show evidence of. glacial action, and 


low interfluve areas with adjacent valleys. Till is 
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restricted to vailey sides. and is rarely exposed from beneath 
lake or fluvial sediments. On the valley sides of this region 
the thin suficial cover. is strongly influenced by local bed- 
rock. This suggests that the surficial material in these areas. 
may be mass movement debris, developed under periglacial con- 
ditions on the hills peripheral to the main valley, rather 
than till. Further research into this problem would be 
necessary to determine whether this material is weathered. 
bedrock, or weathered and/or disturbed till. Evidence of 
intensive periglacial action was found on Tower Hill (Plate. 
4). Almost all of the northwest to southeast oriented valleys. 
of the interfluve areas of the basin are presently occupied 
by bogs, with up to 30 feet of fine clay and gravels. The 
scarcity of sections in till and the poor quality of exposures 
found leaves questions about flow directions and the extent 
of ice cover in the area unanswered. The fluvial environment 
of this part of the study area will be considered in Chapter 4. 
(ii) The second area is bounded by north-south lines 
through McLean Hill and Elbow Bridge (Fig. 20). This area is 
most difficult to describe. The findings about this central 
region of the study area have shown that it is complex. The 
form of the complexities reveal the inadequacy of the data 
collected for a full description of the area both spatially 
and temporally. The basic form of the problem concerns 
boundary zones or limits. One such limit was considered in 
the first part of this chapter, namely the Prairie-foothills 


boundary zone. Others include the western limit of advance 
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of Laurentide ice (Fig. 21), the distinction between two 
lacustrine areas, and the limit of advance of Cordilleran 
ice. 

The area is dominated by northwest to southeast 
ridge and valley forms to the north of Bragg Creek (river) 
and south of the Elbow River. Therefore, this area is con- 
sidered in five units, with the Elbow valley considered in 
Chapter 4. 

(1i-a)The southern area (Region 2A; Fig. 19) of ridge 
and valley is lower in relief than similar topography further 
west. The composition of material covering these slopes is 
Similar to that described on the valley sides in area (1) 
where the true identity of the material was in doubt. 

(1i-b)The ridge and valley topography to the north of 
Bragg Creek town site (Region 2B; Fig. 19) is not as intense 
as that found south of the Elbow River, and tills are more 
clearly identifiable. Valley bottoms are often covered with 
2 to 6 feet of lake sediment with till exposed on spurs and 
flanking ridge slopes. To the north of Bragg Creek town site 
samples of till were taken and originally described as 
Cordilleran in origin. A transect from headwaters to con- 
fluence of Bragg Creek (river) had shown that a small local 
Cordilleran glacier occupied this valley and the tills north 
of Bragg Creek town site were thought to be due to Bragg 
Creek ice and Elbow ice jointly. However, the sediment 
analysis revealed that these tills were not distinctly 


Cordilleran in type. Since the clay-silt fraction was. higher, 
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there is a possibility that the samples involved are from 
Laurentide till. However, throughout. the ridge and valley 
regions and within Bragg Creek valley the influence of bed- 
rock was strong . In this area north of Bragg Creek town 
Site the bedrock is largely Upper Alberta formation which 
consists mostly of soft shale and some sandstone. Thus, 
without further detailed petrographic, mineralogical and 
fabric orientation data, it was not possible to substantiate 
the Laurentide origin. Morgan (1966) and Nichols (1931) both 
described Jumpingpound till in this vicinity as. Cordilleran. 
Nichols (1931) made the suggestion that this till was 
attributable specifically to Jumpingpound Ice. Yet inspection 
of the Jumpingpound basin reveals few ice source areas. 
Therefore, the suppostion by Morgan (1966) that Kananaskis 
and Bow Valley ice spread south, over and beyond Morley Flats 
into the Jumpingpound area seems more feasible. In addition 
researchers from the Calgary base of the Geological Survey of 
Canada (Pers. comm. Dr. N. Rutter) have identified morainic 
forms which are flute-like, and traverse the ridge and 
valley topography of the Jumpingpound basin. These were 
termed flutes but may be ablation moraines attributable to 
melting Bow Valley ice. 

(ii-c)The lacustrine and ground moraines of the area 
peripheral to Bragg Creek town site and extending west and 
south along the Elbow valley define the third area (Region 
2C poe uGeee Lo) s 


All samples from this lake area had greater than 
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98.5% clay-silt. The upper elevation of the lake was deter- 
mined as 3975 feet (the highest elevation at which lacustrine 
samples were found) and the extent of the lake mapped on this 
basis. This lake will hereafter be referred to as Bragg Lake. 
The southern part of Bragg Lake south of the Elbow River is 
dissected by wide, shallow, dry and flat bottomed channels 
which suggests drainage of Bragg Lake to the south east 
through gaps in Sarcee Ridge (Fig. 20). One such gap is 
over 500 feet in total relief and presently carries Priddis 
Creek. The broad, flat lake area bordered by ground moraine 
to the south, Sarcee Ridge to the east, and Last Break Hill 
till area to the north, would seem a feasible pre-glacial 
drainage route for the Elbow River. However, the present 
Elbow River does not take this course but cuts the topo- 
graphic divide between Bragg Creek valley and the Elbow 
valley at a point 1-1/2 miles southwest of Bragg Creek 
town site (Fig. 17). The small northwest to southeast 
oriented valleys east of this point are blocked off on the 
south side by the low moraine traceable to the west for one- 
half of one mile. The Elbow river has breached this moraine 
on the upstream limit of Bragg Creek, whilst an isolated, 
flat, gravel terrace is traceable between the present Elbow 
valley and Bragg Lake. The terraced gravels are weathered, 
covered by Mazama ash and lack good exposures for estimating 
fabric orientation. Imbricate structure shown by pebbles in 
a sand lens, however, indicated flow to the east, directly 
into one of the dry drainage channels of Bragg Lake which led 


to Priddis-Creek. 
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(ii-d)A transect along the forestry road in Bragg 
Creek valley (Region 2D; Fig. 19) revealed a broad, flat 
floored valley with deep morainic infill, and several small 
truncated spurs. This suggested erosion as well as deposition. 
by glacier ice, the source of which must have been in the 
valley headwaters along the eastern.flanks of Moose Mountain 
Dome. End-moraines were traceable down valley to one mile 
east of the forest boundary (Fig. 20). At this point ground 
moraines and lake sediments were encountered as relief became 
low in the lower Bragg Creek basin. Immediately downstream 
of each end-moraine medium to coarse gravels were found, not 
terraced, with no imbrication, and which directly overlay 
bedrock (Plate 5). These gravels were assumed to be outwash 
gravels. North of the Ranger Station Creek (Fig. 20) in 
Bragg Creek valley, are extensive sloughs occupied by lake 
sediments and coarse gravels. The lake sediment uphold active 
bogs. It is thought that these sloughs represent a phase of 
outwash of gravels developed during the final retreat of ice 
from this basin. In the lower Bragg Creek basin, the creek 
has incised deeply into lake sediments and ground moraine 
(Cordilleran) creating widespread and recently active slump 
features (Plate 6). 

(ii-e)Finally, the eastern part (Region 2E; Fig. 19) 
of this central study area, which also forms the western 
limit of area 3 (below and Fig. 19), is a zone not sampled in 
the vicinity of the Sarcee Indian Reserve. It is an area 


dominated by Sarcee Butte to the north and Sarcee Ridge to 
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PLATE 6 


Shows slumping of lacustrine clays at site 131, near 
Bragg Creek. | 


PLATE 7 


Shows banded lacustrine sediments overlying till (lower 
right corner) at Elbow Bridge (Site 1l and 79). 
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the south of the Elbow river. To the east of this zone 
Laurentide tills were identified, whilst to the west 
Cordilleran tills were found. To the east is an extensive 
lacustrine area at a lower elevation than Bragg Lake to the 
west. The degree to which each of the Cordilleran and 
Laurentide Ice masses encroached upon this zone is not 
certain. 

(iii) This area includes all of the study area east 
of the Elbow Bridge. At sites 1l and 79 lake sediments were 
sampled which overlay till in the manner shown in figure 22 
(Plate 7). The lake sediments are of the lower lake, herein 
called Lake Elbow, (although in fact an extension of Glacial 
Lake, Calgary, (Morgan, 1966). The till appeared to be 
Laurentide in origin, as were tills half a mile north (Fig. 
21). These particular till samples fell in the Laurentide- 
Cordilleran. till overlap of. the. factor loading lists. 
Although phase two of the analysis, between the two tills, 
designated the samples from the above sites as Cordilleran 
the initial designation is maintained, because of bedrock 
influences. This is because Bearpaw sandstones are exposed 
locally and are easily eroded, and may therefore have con-. 
tributed coarser material to the tills of this area. Close 
to Hoodoo Ridge in two till sites, granitic pebbles were 
found which further supported the original designation. To 
the east of these sites Laurentide till is exposed on the 
higher elevations above 4000 - 4100 feet (Fig. 21). Below 


this level are lacustrine deposits of Lake Elbow. 
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FIGURE 22 


CLIFF SECTION AT ELBOW BRIDGE 


Shows lake sediments from Lake Elbow overlying Laurentide till. 
The lake sediments are varved ( Plate 17 ) in the top 4 feet, and 


contains buffalo bones. 
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The analysis of the Laurentide tills in phase three 
(components two and four, Fig. 12), suggested two groups of 
this type. The interpretation of the values for each site 
did not, however, produce any understandable relationships 
between the values. This would concur with Morgan (1966), 
who described all Laurentide till in the Elbow basin as 
Erratics Train Ground Moraine. On the other hand the com- 
plexity of till types described by Tharin (1960) cannot be 
envisaged as occurring in this area. In the vicinity of 
Tower Hill to the west of Calgary, evidence of complexities 
of till 7structure was found. Tower Hill, which has long 
smooth convexo-concave slopes on the south side and an exten- 
Sive, flat, very gently dipping eastern summit (of less than 
one degree) (Fig. 20), is everywhere covered by till, except 
on the upper southwest slopes where mass movement has exposed 
bedrock. In addition, in ther benched part of the hill, above 
the lake sediments, periglacial features are common (Plate 8), 
with ball clays (Plate 9) and ice wedges which penetrate bed- 
rock (Plate 10). It cannot be conceived that the bench, cut in 
till and presently dissected by deep misfit river channels, is 
due to a periglacial cycle of erosion such as described by 
Peltier (1950). The till cover precludes pre-Wisconsin forma- 
tion. One suggestion would be the action of Glacial Lake Cal- 
gary. However, this would require the Lake Elbow sediments be 
found up to 4050 feet. In fact lake sediments were only found 
up to 3975 feet, which may indicate mass wasting in the inter- 


fluve areas over 3975 feet. A similar bench is found southeast of 
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Game ice wedge in pebbly till on 
Mae lower Hill at site 24. Gravels 
me underlay the till in the gap at 
the base of the photograph. 


PLATE 9 


Shows what may be a ball clay. The material is till in the 
Same section as above, on Tower Hill. 
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PLATE 10 


Site 51, on Tower Hill. 
Ice wedge penetrates 
bedrock in lower third 
Of photograph. {The 
surficial material is 
faiede. 


PLATE 11 


Shows bench in the centre horizon, with two stands of trees 
on it. Location is on a S.E. extending ridge from Tower 


Ridge. 
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Tower Ridge (Plate 11) projecting into the lake region at the 
4000 foot level. To the east and west of Tower. Hill are over- 
flow channels at 3725 feet and 3870 feet respectively. In 
each case deltas are built up on the south side indicating 
overspill from the Bow valley into Lake Elbow at two different 
time periods. To the east of Tower Hill and interbedded with 
the delta on the east side, is till which forms a notable 
north to south ridge. This is thought to be an end moraine. 
The land elevation to the south of the Elbow in this part of 
the study area is generally below the 4000 foot level. 

Conclusion for the Post Tertiary Geomorphic 
History of the Elbow Basin 

The conclusion pertaining to the possible sequence of 
events in this area are now presented. 

Preliminary investigations suggest Laurentide ice 
advanced as far west as Sarcee Butte and Sarcee Ridge. At 
this time, Bragg Lake was formed as the Priddis Creek gap and 
the present Elbow Valley were blocked off. Subsequent retreat 
of Laurentide ice in a north-east direction allowed Bragg Lake 
to drain south-east through Priddis Creek before increased 
volumes of meltwater eventually cut the present course. The 
Laurentide ice front must have been extensive, and when it 
was in retreat meltwaters formed in increasingly extensive 
ice front lake - Glacial Lake Calgary. Subsequently, 
Laurentide ice either re-advanced or stagnated on the eastern 
flanks of Tower Hill. The lake level rose to approximately 


4050 feet when the two benches (Fig. 20) were cut in ground 
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moraine. With continued retreat the lake level dropped as 
water drained either directly south or through the present 
Elbow valley south of Tower Hill. Exposure of the benched 
areas in a lake environment created intense periglacial 
action. This lake sediment, deposited on high interfluve 
areas, was washed off and till exposed as water drained away 
(Fig. 21). This would explain why lake sediments are not 
found up to 4050 feet. As the Elbow - Bow valley watershed 
was exposed overspill from the Bow cut the higher spillway 
(Fig. 23). The continued Stagnation of Laurentide ice and a 
drop in lake levels created Lake Elbow, until eventually water 
drained south again on the east flank of Tower Hill. A minor 
advance is illustrated in Figure 24, where till overlays 
disturbed gravels and is overlain in turn by very coarse and 
massive bedded deltaic gravels (Plate 12 and 13). Subse- 
quently, the lakes were completely drained as ice retreated 
again. A paleo-soil found close to the delta (Plates 14 and 
15) could not be interpreted without further information from 
within the vicinity. However, the following additional points 
are relevant to the conclusion. 

(1) The lake Elbow sediments contain many pebbles 
which were probably ice rafted into the area. 

(ii) Buffalo bones were found in the top 18 to 24 
inches of the lake sediments, indicating a prolonged phase of 
existence of Lake Elbow as a distinct lake when 'life' could 
return to the area. In addition, it would place the final 


retreat of the ice within very late Wisconsin times. 
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FIGURE 23 


Silt and fine sand. No imbrication. 


SKETCH FROM SITE 36 


On east flanking slope of spillway near delta to the 


west of Tower Hill 


The section is close to the floor of the spillway. The distortions in the till and in 


the coarse gravels are small jateral scour channels. 
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PLATE 12 


These gravels are deltaic. 
There is no visible imbri- 
cation. The sides of the 
Square are each one foot 

in length. The section is 
found at McPherson's Quarry, 
just west of site 25, in 
S.W. Calgary. 


PLATE 13 


Shows till enclosed in gravels which form part of the delta 
at site 25. The road is Sarcee Road in S.W. Calgary. 
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PLATE 14 


Paleo-soil buried by mass waste material near site 25. 


PLATE 15 


Continuation of the paleo-soil shown in plate 14, with 
lacustrine sediments overlying gravels. 
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(iii) The lack of sediments above 3900 feet is due to 
two factors. The first is that the interfluve areas are 
barely above this elevation, and hence deposition of material 
in these areas would be limited. Secondly, in the high 
relief areas around Tower Hill mass movement is evident. The 
proposed type of movement is by mudflow, as a section taken 
at 3874 feet at site 41 (Fig. 25 and Plate 16) shows overlap- 
ping shallow channels on a 15 degree slope in very fine clay 
material. This would also explain local, very stony patches 
of till in which the clay-silt content was low (around 30% to 
40%) for Laurentide till. Mass movement by soil creep is held 
responsible for exposing till around Tower Ridge. 

(iv) The form of the bench on Tower Hill suggests 
very shallow water. This is the only place a shoreline may 
be located. Elsewhere, the drop in lake level must have been 
fairly rapid. The misfit valleys of Tower Hill may also be 
accounted for by lake levels which covered the bench on which 
they are found. 

The chronology pertaining to the Cordilleran ice was 
indeterminable on the basis of data collected. The relation- 
ships of the post-glacial environment in the Elbow River will 


be discussed in Chapter four. 


gor ee PU iga N ea a 
ot eu ef test 00E evods ataomiber 30 oat of? ( 
ais aseis eveltretal end tens at textt edt 87 
{sizetam 20 moksisogeb eoaedt ‘bow inks sveke eta : . eh | 
pia sate mk se bacnae bad bmbs 24 bluow ‘asoxs 8 
edit .toebive ef $cemovont geem LIth s9Wor babors esors | 
netet nolsose s as <wolthum yd ak Srnemevom: 0 equ? 6 
-gslisvo awete (el esale Bas es. ptt) - all 36 toot wees 
- ypfo oni yxev ai egole ob-ipeb @f 5 mo abeniedo eesti 
aerotsq yore YxeV ,Leoot nisiaxs oals hivow aidt sited: on 
o+ #0€ Brwors) wold @bw JSaezn00 tiie-yeis ods dotdw at sa Mir 
bled at aqesxvo Lhom yd tnemeyor eaam mage; aap ai | 
-epbin rzswoT bewets Iles ei ieoute 10% : 
eteeppia ILih sewolT ao domed ent To iy ont sii! aS ‘i 


% rm 


ysm onifexode s sosig ybito git ‘el vid? .19Few 
maed evssi eum fsvel sxsi mt gorb edt voxodwen tt vaneail 
sd oels yaa Sidi 1awor to ayeliav titeim anit eauend ‘ ates : 
dekiw-no doned ett bexsveo dotnw alevel ann i ees. 6 
| Nae ous 

ssw opi netell bexoS Sas oF priaisitee wiltondges edt 
~noiiate: sd .Bedoetfes agab Zo ehesd edt no eldentmsosabat 
fLiw 1avid wodld edd me sscunaaimion Rptosty sant ont to weiss 


\ saad 2 > a ome ed 


109 


Ground Surface 


Ice Wedge 


4.25 feet 


FIGUREMZS 


SKETCH OF SECTION SHOWN IN PLATE 16 


The section is cut on a concave lower slope on the S.W. side of Tower Hill. 


A Silt and clay 

B Dark grey clay 

Cc Silt 

D Sand and clay 

E Silt and clay 

F Dark clay 

G Silt, clay and sand 

H Fine sand, silt and clay 

I Dark clay with pebbles up to 8mm 
J Silts, sands and gravels . 
K Clay 
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PLATE 16 


Site 40. Mudflow patterns 
are partly discernable in 
a clay rich material 
(probably lacustrine) 
containing thin bands of 
small pebbles. 


PLATE 17 


Photograph of fossils collected from lacustrine sediments at 


Site 6, cat .Connop. Creek 
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CHAPTER IV 
THE ELBOW RIVER AND VALLEY 


The behaviour of the Elbow River within the present 
valley is discussed in this chapter. The time period 
studied is from late Wisconsin times to the present. 

The upstream section of the Elbow valley is 
dominated by pedence constrictions between which are 
numerous unpaired terraces. This was established from 
approximate altimetric measurements of all the major 
terraced surfaces upstream of Bragg Creek, taken from 
aerial photographs at the 1:2640 scale. South of the 
river and above the highest terraces the land is relatively 
flat where ground moraines border the fluvial debris. [In 
the downstream section the valley is cut in lake sediments. 
Bedrock is rarely exposed, and no terraces are found. 
Occasionally the river has cut benches in the lake sediments. 
The point of change between the two sections is near Bragg 


Creek town site. 


The Elbow River Upstream of Bragg Creek 
Upstream from the road bridge at Bragg Creek the 
behaviour of the Elbow River is subject to bedrock control 
in two ways. The first concerns constriction of the valley 
where bedrock outcrops in the valley side, thereby restrict- 


ing the continuity of the river pattern within the floodplain. 
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The second bedrock control concerns the exposure of bedrock 
in the side of the channel. Associated with each bedrock 
control is a local steepening of the channel gradient, 
except at the bedrock outcrop immediately upstream of Bragg 
Creek (at the Provincial Park). In most cases, the steep 
gradient is found in the region of the outcrop or on the 
upstream side (Fig. 26). This eee suggest progressive 
headward retreat and the removal of the influence of each 
local bedrock base level. At the same time, the majority 
of these points continue to influence the slope of the river. 
The exception is found upstream of Bragg Creek at the 
Provincial Park (Fig. 2) constriction. In this case, soft 
carbonaceous shales of the Upper Alberta formation have been 
more rapidly eroded, such that adjustment of slope to the 
downstream local base level at the road bridge in Bragg 
Creek is now taking place. The bedrock is very close to 
the surface in the Provincial Park constriction, but ene 
river has aggraded sufficiently to produce a braided 
channel. The aggradation is thought to be quite recent. 

The reaches between each constriction upstream of 
the Provincial Park are incised into thick deposits of 
coarse fluvial material. In some cases lacustrine and/or 
till deposits are exposed in the valley side beneath the 
fluvial material (Fig. 26). One such occurrence is shown 
in Figure27, at Connop Creek. In this case a collection of 
small fossils were found in clay rich lake sediments, 


consisting of gastropods and occasionally brachiopods 
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FIGURE 27 


ROAD CUT SECTION AT CONNOP CREEK 


Shows coarse fluvial gravels (A) overlying lake sediments (B and C). 
The till (D) is Cordilleran. The lake sediments ore part of a small local 
lake ond are varved. The thin bands (C) are black, clay rich and contain 


small fossils (Plate 17) 
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(Plate 17). Another site at which similar fossils were 
collected was found at the top of the gorge 300 yards west 
of Canyon Creek. In this case lake sediments largely 
composed of silt and fine sand overlay till and bedrock. 

At another site half a mile east of the Kinsman Camp 
constriction, at site 83 (Fig. 2), evidence of at least 

one re-advance of the Elbow valley ice was found (Fig. 28). 
The section is taken from a 50 foot river cliff on the south 
Side of the floodplain in which erosion at the cliff base 
has caused slumping of the cliff. The till everlain on 
the west side by fluvial gravels, and overlaying gravels on 
the east side, contain a large number of well rounded 
pebbles (Plate 18). This suggests that the buried and 
disturbed gravels represent a phase of outwash prior to 
re-advance of the ice. 

Despite the above localized complexities of the 
stratigraphy of the upstream section of the Elbow River, 
there is one unifying factor. The coarse gravels associated 
with the highest terraces have buried or replaced most of 
the earlier deposits. The same gravels are characterized 
by the Lack of Stratification. Sand lenses = (Plate 20) are 
common, and a few pebble orientation fabrics taken vary in 
direction and strength of the major orientation. The 
uniformity of the material found within the terraces is 
striking (Plate 19). Grain sizes range in excess of 256mm, 
and boulders up to 18 inches are not uncommon. As shown in 


Fig. 9, the percentage finer than 0.063mm is small. The 
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Slumped 


W ater surface of Elbow River 


FIGURE 28 . 


SKETCH OF SECTION IN RIVER CLIFF 


AT SITE 83 


Shows till G overlying disturbed gravels A with gravel B above 
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PLATE 18 


Shows gravels 'A' (Fig. 28) in the lower left. Till contents 
include a mixture of well rounded pebbles. (Till = Till I 
figure 28). 
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PLATE 19 


Section from garbage pit 
at site 3, near the Elbow 
Ranger Station. Shows 
medium coarse gravels 
poorly bedded. 


PLATE 20 


Shows sand lens in gravels exposed at site 3. Current 
flowed from left to right. 
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19 
grain size distributions fall into two groups. The first 
group is unimodal with a median grain size between 16mm and 
32mm. The second group is bimodal with peaks between 4mm 
to 8mm, and 32mm to 64mm. The bulk field samples referred 
to in Chapter 2 were taken near Elbow Ranger Station 
(site 3) and half a mile east of the Kinsman Camp (site 83). 
The results of comparing samples taken from the horizontal 
and the vertical in each section showed negligible 
differences in the total percentage contribution of 


material greater than l6mm, as the figures below illustrate; 


Site Horizontal Vertical 
3} Silas 6 5229 % 
83 5Ba3. Ss Borea2 - & 


The samples from the horizontal (Plates 21, 22, and 23) might 
be considered as representing distributions of material in 
space at one point (or one depositional phase) in time. 

The samples from the vertical (Plate 23) might represent the 
mean depositional characteristics in time, either in one 
phase or between several phases of deposition. In essence 
therefore, the above suggests that the depositional environ- 
ment was consistent in time and space. The greatest thickness 
of gravels is approximately 100 feet, found locally between 
the Elbow Ranger Station and Beaver Flats. Throughout. the 
upstream section the valley sides are benched. The fluvial 
material beneath all the benches is similar to that described 


above. In addition, the benches are cut in till and bedrock 
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PLATE 21 

Ee 

Bee) Grid overlying section on 

Pe the horizontal--cleared to 

ae Way a depth of 9 inches. Site 
‘Meets. Sampe 400. 1970. 


PLATE 22 


As in plate 21, 
showing coarser 
gravels. Site 

83. Sampe 300. 
LOO? 


é 
: | 
ae 
s | PEATE =25 
Taken at the same site as 
¢ 
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<8 lack of fines (less than 
8mm) . 


plate 22, but in a vertical 
section. Note the apparent 
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at different locations and varying elevations. The 
conclusion is that the gravels accumulated rapidly, with 
each successive retreat of the Elbow ice from constriction 
to constriction. As Figure 29 indicates, the greatest phase 
of outwash occurred when ice had retreated as far as Beaver 
Flats. 

In a few sections a form of bedding structure is 
obvious. The bedding is distinguishable in terms of the 
matrix for pebbles greater than 8mm. The distinctiveness 
of the bedding is therefore in the clay, silt and sand 
fractions. One such case is found in Bragg Creek just below 
an end moraine at site 106 (Plate 24). In the higher 
terraces upstream of McLean Creek (Fig. 29), the bedding is 
less obvious, though at site 135 a more strongly bedded 
section was found. This contained one bed (Plates 25 and 
26) which had no material less than 16mm in size. This may 
indicate a form of seasonal persistence of flow with seasonal 
high flows of considerable discharge. The changing nature 
of the matrix suggests extensive channel migration ina 
braiding system. In all, this would indicate the initial 
aggradation of the valley infill as an outwash plain or 
sandur (Plates 27 and 28). This is because valley reaches 
between successive constrictions seem to have individual 
terrace and bench sequences. The only exception is found 
near the Elbow Ranger Station where the highest terrace 
(number 12, in Fig. 29) at this point of the valley extends 


downstream across the Elbow Ranger Station constriction. 
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Silt and sand matrix 


Cobbles and boulders. 


Gravel and sand matrix. 


Cobbles, 


PLATE 24 


Site 106. The ruler transgresses two types of matrix in the 
gravels. The finer matrix is associated with re-excavation 
or re-exposure of till for the transportation of fines. 
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PLATE 25 


General view of section at site 135. 


PLATE 26 


Shows gravel bed at the 
base of the ruler in 

which there is no material 
less than 16mm in size. 


Arrows on plate 25 locate 
the position of the ruler 
in site 135, in which this 
photograph was taken. 
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PLATE 27 


View east from position 200 yards west of site 149 (west of 
Canyon Creek). Note extensive gravel terraces in the base of 
the valley. 


PLATE 28 


View west from site 149. Shows Prairie Mt. in the left centre, 
and the upstream end of the extensive terraces shown in plate 
Pe es 
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Upstream of this point till is found in isolated sections 
below the terrace gravels. The top of the till near the 
Elbow Ranger Station is almost at the same level as the 
bedrock in the constriction at this point in the valley. 
This may indicate an initial phase of degradation in the 
glacial infill prior to accumulation of outwash. Hence, 
the sandur may have developed on a plain which was graded 
by the river between the bedrock knickpoint at the Elbow 
Ranger Station, and the snout of the glacier near Beaver 
Flats. At the point in time when this major aggradational 
phase terminated, the assymetry of the Elbow River and the 
bedrock valley was most probably established. This is 
because the river became independent of the control of a 
bedrock valley. 

Channel incision occurred subsequent to the climax 
of aggradation of the valley. As shown in Figure 29, the 
degradation was fairly continuous, as suggested by the 
unpaired nature of the benches. The problem of correlating 
benches is heightened by a relatively steep gradient in most 
cases upstream of the Kinsman Camp. Hence, evaluation of any 
notable periods of incision is not possible at this time. 

The deep gorge traceable for one mile below Beaver 
Flats is cut through gravels and then into bedrock. When 
the gorge was cut is not directly determinable. If the 
gorge is pre-glacial, it seems improbable that the same gorge 
was re-occupied without considerable lateral wandering of 
the river. The latter would not permit the side of the 


gorge to extend vertically into the gravels. However, 
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within the gorge active rock falls are found and the Elbow 
Falls shows signs of retreat upstream. For example, scour 
marks are found in front of the falls. This suggests post- 
glacial formation of the gorge. The only dating factor for 
the change from an aggradational to a degradational fluvial 
environment is that Mazama ash is found on the highest 


terraces upstream of the Elbow Ranger Station. 


Summary 

At the present, upstream of Bragg Creek town site, 
the river appears to be adjusting to local changes in base 
level. Hence, between Bragg Creek and the Kinsman Camp, 
degradation is active as the river adjusts to the lower 
base level at Bragg Creek road bridge. Upstream of the 
Kinsman Camp, the local base levels still operate to some 
extent, although it is felt that a point is being reached 
where long term flow characteristics will determine the 
aggradation and degradation response of the river. An 
additional factor will be changes in the courses of the 
river. At present the river is located in mid-floodplain 
or erodes laterally into a bedrock valley side. Should 
the river course adjust to erode into the benched gravels 
then more active aggradation may occur. In essence therefore, 


a quasi-equilibrium may have been attained, with a greater 


tendency for degradation. 
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The Elbow River Downstream of Bragg Creek 

In the downstream section below Bragg Creek road 
bridge, the major controls are the valley side material, the 
bedload material, and differential transportation of material. 
The initial establishment of the valley postdates Bragg 
Lake and Lake Elbow. The depth of coarse gravels (which 
are as coarse as upstream of Bragg Creek) is in excess of 
30 feet in places above bedrock (Hollingshead, 1968). It 
is probable that the valley was established synchronous 
with the building of the outwash plains. Subsequent 
transport of coarse gravels downstream increased the depth 
of the valley by concentrating flow with a heavy bedload. 
In theory therefore, the influx of ccarse debris caused 
aggradation once bedrock had been reached. The degree of 
control of base level at the confluence with the Bow River 
is difficult to establish. The Bow River has incised 
considerably, but Stalker (1961) states that the modern Bow 
River is higher than the pre-glacial Bow River, which 
suggests continuing degradation at the Elbow's confluence 
with the Bow. Evidence of initial adjustment of the Elbow 
River to a lower base level is found in the steep valley 
sides directly flanking the present floodplain. Yet, as 
mentioned above, 30 feet of coarse gravels have accumulated 
in this valley. It appears that change in the bedload 
material is the major factor to be considered in the 
evaluation of the aggradational-degradational behaviour for 


the immediate post-glacial period. Hence, the following 
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concerns the recent behaviour of the Elbow River. 

The floodplain is widest near the Elbow Bridge. 

At present the river in this stretch and for most of its 
course below Bragg Creek is on the north side of the 
floodplain. However, sixty years ago the main stream was 
on the south side in this stretch, according to personal 
communication with the local inhabitants. Also within the 
Same area, aerial photographs show that white spruce and 
lodgepole pine stands differentiate small differences in 
elevation of channel bars. From field observations these 
differences reflect increased incision of the river, as 
well as a tendency for migration of the less intensely 
braided open channel areas (Plate 29). The channel pattern 
is essentially an irregular meander, which during June 
high flow periods becomes markedly braided. The remnants 
of the braided channels fall in zones as mentioned above. 
A noticeable occurrence is the association of cross valley 
benches with straight, less braided channel reaches. This 
suggests local aggradation and degradation reaches. 

In 1968, Hollingshead noticed a marked decrease in 
median grain size close to Twin Bridges, which continued 
downstream. It is unlikely that this change is related to 
the Glenmore Reservoir, which is now the base level of the 
Elbow River. The cause of the change in bedload 
is probably attributable to flow and selective transportation, 
rather than the reservoir. Since the finer bedload extends 


into the reservoir, the conditions causing the formation 
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PLATE 29 


View §.W. from site 79 shows the floodplain of the Elbow 
River upstream of Elbow Bridge. Note the vegetation on the 


floodplain. 
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of this reach must still operate. Hence, this accumulation 
of finer gravels is continuing. The lower reaches of the 
downstream section are steeper and more irregular than 
might be anticipated. The local steepening of the gradient 
approximates the gradient at the cross valley benches 
mentioned earlier. A possible cause of the gradient changes 
may be limits of advance of debris loads attributable to 


phases of high flow in neo-glacial times. 


Summary 

It would seem that below Bragg Creek the coarse 
gravel reaches are being degraded as finer gravels are 
carried downstream, and clay and silt accumulate in the 
Etoodplain. “In the lower’ reaches aggradation is occurring. 

The behaviour of the Elbow River within the study 
area is changing as present flow conditions cause increasing 
selective transportation. As a result of this action and 
bedrock controls upstream of Bragg Creek, changes in quasi- 
equilibrium reaches may contribute to changes in strongly 


aggrading reaches. 
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CONCLUSION 


This study investigates morphological and sediment 
characteristics, and analyses the behaviour of the present 
Elbow River and its basin. The major objective was to 
provide a basis for assimilation of the earlier, more 
specific studies, and to indicate useful further research. 

The information upon which this study is based 
falls into three categories. The first concerns sediment 
samples which were collected and analysed for the grain size 
content. The results of this analysis were then used to 
determine the reliability of the original diagnosis of the 
sediment type of each sample made in the field. Principal 
component analysis was used to classify sediments according 
teegrainsize scontent 4° Justificat sonvforrthe 4uselCof@principal 
component analysis in this context is presented in Chapter 
2. The second category concerns the collection of 
morphological information in the field, from aerial photo- 
graphs, and from topographic maps. The third category 
involves the collation of information, from previous 
studies, concerning the geologic history of the study area. 

The findings of this study may be considered in 
four parts. The first is pre-Pleistocene geologic history, 
the second Pleistocene gromorphology, the third part 


behaviour of the Elbow River, and fourth the granulometric 


analysis. 
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The analysis of summit peak data from within the 
Elbow basin indicated that the topography of the foothills 
is closely related to geologic structure, and that the 
topography of the Prairies is more independent of structure. 
In addition, it has been shown that a sub-aerial environment 
existed in the study area prior to earth movements of the 
Laramide orogeny, and this environment persisted throughout 
the Tertiary era. It has proved difficult to conclude 
whether drainage is antecedent to geologic structure or 
superimposed from a mid-Tertiary peneplain in the foothills. 
The drainage of the Prairies, marginal to the foothills, 
is more closely related to peneplanation forms of a 
Pertiary=age. 

The Pleistocene geomorphology of the basin is 
attributed to the effects of the Wisconsin glaciation. The 
two ice masses which affected the study area are Cordilleran 
and Laurentide in origin. The extent to which each ice mass 
covered the area is inferred from two separate parts of 
the study area. The eastern part is associated with 
Laurentide ice deposition in the prairies, while the western 
area is associated with erosion and deposition in the 
foothills by Cordilleran ice. The dividing line between 
these two parts of the study area is associated with the 
western limit of advance of Laurentide ice, which is 
thought to be along a line through Sarcee Butte and Sarcee 
Ridge. The formation of two ice front lakes is attributed 


to (a) the limit of advance of Laurentide ice (Bragg Lake), 
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and (b) to a minor advance stage during the final retreat 
of this ice (Lake Elbow). In the foothills the glacial 
deposits are attributed to Cordilleran ice which originated 
within the Elbow basin, upstream of Beaver Flats. A smaller 
glacier is thought to have occupied Bragg Creek valley with 
a source area on the eastern flank of Moose Mountain Dome. 
The main depositional feature of the foothills in the study 
area is a thick infill of fluvial sediments in the Elbow 
River valley, upstream of Bragg Creek. These deposits are 
thought to be part of an outwash plain. 

The present Elbow River, downstream of Bragg Creek, 
appears to be degrading in the coarse gravels, as finer 
gravels are carried downstream. In the lower reaches, above 
the Glenmore Reservoir, aggradation of the finer gravels 
occurs. Upstream of Bragg Creek selective transport of 
sediment, and the frequent occurrence of bedrock control, 
cause local quasi-equilibrium states of the longitudinal 
profile. As bedrock control is minimised, degradation 
extends upstream until a new quasi-equilibrium state is 
reached. This is believed to account for the behaviour of 
the Elbow River at the Provincial Park, near Bragg Creek. 

The analysis of grain size data has shown that 
principal component analysis can be used to discriminate 
between sediment types. In addition, this statistical 
technique identified the obscure cases where additional 
geological information was required to substantiate the 


original sediment type allocation. 
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Future Research 

In the author's opinion the determination of the 
geomorphic history of the Elbow basin depends upon more 
detailed examination of the surficial deposits. The 
systematic examination of morainic material, deposits on 
hillslopes in the foothills, lake sediments, and the 
relationship of ash deposits to terrace levels are needed. 
The study area, which includes part of three major 
topographic regions in wleeHern Alberta, may yield valuable 
information concerning the pre-Wisconsin glacial events, 
limits of advance of Laurentide ice, and the modification 
of pre-glacial drainage. 

A systematic investigation of the Elbow River is 
warranted in view of the change in sediment regime upstream. 
The primary features requiring investigation are sediment 
transport rates, hydrological characteristics, channel 
behaviour, and river patterns. 

Principal component analysis is a useful tool in the 
field of environment discrimination in the spatial context. 
It is felt that the procedure outlined in Chapter 2 could 
be applied to a smaller area, in which a more controlled 
sample plan for the collection of sediment samples could be 
employed. This would permit a more reliable understanding 


of the statistical and spatial aspects of this analysis to 


be achieved. 
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APPENDIX 1 


PROGRAM LISTING FOR CALCULATION AND PLOTTING OF 


GRAIN SIZE DATA 


This program was prepared by G. C. Seagel, in 


December, 1970. The program calculates the percentages, 


cumulative percentage and an expression for the logarithm 


of the cumulative percentage. The latter is generated from 


natural logarithms to enable the plotting of the curve. 


The program is easily arranged so that Folk and Ward (1957) 


descriptive statistics may be generated. 


SIEVE me 


PHI = 


The program 


mesh size in millemetres. 

phi equivalent for sieves used. 
identification number of each sample 

( neta Sbeormac.) 

site number ( in 'A' format). 

weight of material retained by sieve 
X, and associated with sieve X + 1, as 


sieved finer than sieve size X + l. 


is presently arranged to take and 


compute for 19 consecutive grain size classes. 
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The Program: 


DIMENSION SIEVE (20) ,PHI (20) ,WT (20) ,CUMWT (20), 
PCWT (20) ,CUMPC (20) ,XL 


1(20) ,XLOG (20) ,XNORM (20) ,MIN (20) ,MAX (20) ,NEWIN (20), 
THIST (20,101) ,LP 


200420 , 101) 7 GRP (7) 
INTEGER EYE,CUPLT (20,101) 
INTEGER, BL/*!/,FD/'+'/ 
DOM 15 29 
READ (5,1001) SIEVE(J) ,PHI (J) 
1001 FORMAT (F7.3,F4.1) 
1 CONTINUE 
99 READ(5,1002) NSAMP,LSITE 
1002 FORMAT (A4,A4) 
READ (5,1003) (WT (J) ,J=1,19) 
1003 FORMAT (10F8.2/9F8.2) 
po 100 I=1,19 
DO LOL J=1,101 
IHIST (I,J) =BL 
CUPLT (I,J) =BL 
101 LPLT(I,J)=BL 
100 CONTINUE 
G CALCULATE CUMULATIVE PERCENT 
TOTWT=0.0 
DO .10 J=1,29 
TOTWI=TOTWT+WT (J) 
10 CUMWT (J)=TOTWT 


CUM2=0.0 


ant 


, (08) HDD, (0s) mm, (08) f er 


, (os) wIWaM, (OS) XAM, (OS) WIM, (ony 


G1, (1QL, os)varur nen oe 


ra “ey a0, (iota) 
Pebiesieipo ge ncire pois ale 


i at) BY n 


Had | co atin penny 
(o) ee 
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np ‘PAGGRS AVERAIUMUD STALUDIAD 
amas intel nity Saal 


e i) 4 
to, vhs 


iy cohen ied: sp Haha 
| if 


Et 


ee) 


22 


23 


DONPIAJ=1 ,19 

CUM= (WT (J) /TOTWT) *100.0 

PCWT (J) =CUM 

CUM2=CUM2+CUM 

CUMPC (J) =CUM2 

CALCULATE GROUP PERCENTS FOR SEVEN GROUPS 
DOM SUIS 1,77 

GRP (I)=0.0 

GRP (1) =PCWT (1) 

GRP (2) =PCWT (2) +PCWT (3) +PCWT (4) +PCWT (5) 
GRP (3) =PCWT (6) +PCWT (7) +PCWT (8) +PCWT (9) 
GRP (4) =PCWT (10) +PCWT (11) +PCWT (12) +PCWT (13) 
GRP (5) =PCWT (14) +PCWT (15) 

GRP (6) =PCWT (16) +PCWT (17) 

GRP (7) =PCWT (18)+PCWT (19) 

CALCULATE LOG. CUMULATIVE PERCENT 

DO 20 J=1,19 

XL (J) =CUMPC (J) *10000.0 

XLOG (J) =ALOG (XL (J) ) 

XNORM (J) =( (XLOG(J) *10.0) -130.0) *10.0 
IF (XNORM (J) .LT.0.0) XNORM (J) =0.0 
LOAD=IFIX (XNORM (J) ) 

REG=XNORM (J) -LOAD 

IF (REG-0.5) 22,23,23 

MIN (J) =LOAD 

Go ro 20 


MIN (J) =LOAD+1 
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CONTINUE 
LIST SIEVE SIZE, PHI, PERCENT WT, CUMULATIVE PERCENT 


WIr,LOG. OF CUMULATIVE PERCENT, AND SEVEN GROUP 
PERCENTS. 


WRITE (6,2002) NSAMP, LSITE 


FORMAT ('1'////,40X,'GRAIN SIZE ANALYSIS'//10X,' 
SAMPLE NUMBER = ',A 


14,10X,'SITE NUMBER = ',A4//10X,'SIEVE SIZE IN M 


Me 5x, PHL: VALUES, 


210X,'"PERCENT WT',13X,'CUMULATIVE PERCENT! ,5X,'L 


OG. OF CUMULATIVE P 


3ERCENT' ) 


DORSOTI=1,19 


WRITE (6,2003) SIEVE(J) ,PHI(J) ,PCWT (J) ,CUMPC (J), 
XNORM (J) 


FORMAT (/10X,F7.3,15X,F4.1,15X,F8.4,15X,F8.4, 
15X,F8.4) 


CONTINUE 
WRITE (6,2050) (GRP(I) ,I=1,7) 


FORMAT ('1'////10X,'GROUPED DATA BY PERCENT WT. 
, AFTER WENTWORTH'// 


1/10X,'SILT-CLAY' ,11X,F8.4//10X, 'FINE SAND!',11X, 
F8.4//10X,'MEDIUM S 


2AND' ,9X,F8.4//10X,'COARSE SAND',9X,F8.4//10X,'G 
RAVEL',14X,F8.4//10 


3X, 'COBBLES' ,13X,F8.4//10X,'BOULDERS' ,12X,F8.4) 
PLOT HISTOGRAM 
CONVERT TO WHOLE NUMBERS 
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KOLD=IFIX (PCWT (J) ) 
REK=PCWT (J) -KOLD 
IF (REK=0.5)41,42,42 
41 MAX (J)=KOLD 
GO TO 40 
42 MAX(J)=KOLD+1 
40 CONTINUE 
DO 51 J=1,19 
K=MAX (J) 
51 IHIST(J,K)=FD 
WRITE (6,2010) NSAMP,LSITE 


2010 FORMAT ('1'///20X,'HISTOGRAM OF PERCENT WT. I 
N EACH GRAIN SIZE CLAS 


1S'//15X,'SAMPLE NUMBER = ',A4,10X,'SITE NUMBER = 
' A4///56X,'PERCE 


2NT WEIGHT') 
WRITE (6,2011) 


2011 FORMAT (//17xX,'0",8X,"'10',8X, "20" , 8XVEspeyexkr 
"40" 7 SoG, 50°79, 8X7. 460" 45 


1X70 CSK, 280! , 8x v9 0 ea COL dre ae POLLS Reena yy 
DO 54 I=1,13 
WRITE (6,2012) SIEVE(I), (IHIST (I,K) ,K=1,101) 

2012 FORMAT(14X%,! !'/4% 82 .3725)4- 4 ron LOLA) 

54 CONTINUE 

DO 59 I=14,19 
WRITE (6,2912) SIEVE(I), (IHIST (I,K) ,K=1,101) 

99012 FORMAT (14x%,° '/145,° 1/74X,87-3,25:4- arasnyls)) 


59 CONTINUE 
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PLOT CUMULATIVE PERCENT CURVE. 
CONVERT TO WHOLE NUMBERS. 
DO-60.J=1,19 
HOLD=IFIX (CUMPC (J) ) 
RES=CUMPC (J) -HOLD 
IF (RES-0.5)61,62,62 
NEWIN (J) =HOLD 
GO To 60 
NEWIN (J) +HOLD+1 
CONTINUE 
DO: 65 J=1,19 
K=NEWIN (J) 
CUPLT (J,K)=FD 
WRITE (6,2020) NSAMP,LSITE 


FORMAT ('1'///20X,'PLOT OF CUMULATIVE PERCENT CURVE' 
//15X,'SAMPLE N 


1UMBER = ',A4,10X,'SITE NUMBER = ',A4///56X, 'CUMULAT 
IVE PERCENT WET 


2GHT' ) 


WRITE (6,2011) 

DO 68 I=1,13 

WRITE (6,2021) SIEVE(I) , (CUPLT (I,K) ,K=1,101) 
FORMAT. (LAX, 1/ 4x, 67.3, 2X, = =, 3X, 1014) 
CONTINUE 

pO 69 I=14,19 

WRITE (6,2921) SIEVE(I) ,CUPLT (I,K) ,K=1,101) 


FORMAT (14X," "JAK? “74K B.S peks = “ak, LOTAL) 
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69 CONTINUE 
PLOT LOG. CUMULATIVE PERCENT CURVE. 
DOKTE Teh, LO 
K=MIN (J) 
71 LPLT(J,K)=FD 
WRITE (6,2031) NSAMP,LSITE 


2031 FORMAT('1'///20X,'LOG. PLOT OF CUMULATIVE PERCENT 
T CURVE'//15X,'SAM 


1PLE NUMBER = ',A4,10X,'SITE NUMBER = ',A4///56X, 
"LOG. SCALE CUMULA 


2TIVE PERCENT' ) 
WRITE (66,2032) 


OSD KORMAT(//L7X, 50%, 78X, 1004 11 0G. VAEUE 707x, 
' EO OE n) 


DO, 79 1=1713 
WRITE (6,2033) SIEVE(I) ,(LPLT (I,K) ,K=1,82) ,XNORM (I) 
2033 FORMAT(14X,' '/4X,F7.3,2X,'- ',3X,82A1,10X,F8.4) 
78 CONTINUE 
DO 79 I=14,19 
WRITE (6,2933) SIEVE(I), (LPLT(I,K),K=1,82), XNORM(T) 


5933 FORMAT (14X,' '/14X," '/4X,F7.3,2X,'— %,2%782A1,10X, 
F8.4) 


79 CONTINUE 
PUNCHED OUTPUT 
WRITE (7,3001) NSAMP, (GRP (I) ,1=1,7) 


3001 FORMAT (A4,7F8.4) 


WRITE (7,3002) NSAMP , (PCWT (J) ,J=1,19) 
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LEZ 
FORMAT (A4,9F8.4/10F8.4) 
WRITE (7,3002) NSAMP, (CUMPC (J) ,J=1,19) 
CHECK FOR ANOTHER SAMPLE TO BE PROCESSED. 
READ (5,1005) EYE 
FORMAT (I1) 
TEP .Gl .0)GO.,rO 99 
STOP 
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APPENDIX 2 


PROGRAM LISTING FOR RANK ORDERING OF GROUPED DATA 


BY THE FIRST VARIABLE 


This program was prepared by G. C. Seagel, in 
December, 1970. The program was prepared specifically for 
use with granulometric data. The format statements are 
Geciiy, Changed to accommodate a Greater range Of data. 

N is the number of vectors of observations. 

WI is the vector of observations of length J (in 

this case J = 8). 
S is dimensioned equal to J. 


SWT is dimensioned equal to I equal to N. 


The program: 
DIMENSION WT (300,8) ,S(8) ,SWT (300) 
WRITE (6,200) 


200 FORMAT('1'///46X,'RANK ORDER OF GROUPED DATA & 
PERCENT WEIGHT'//) 


READ (5,100) N 
101 FORMAT (13) 
READ (5,101) (WI(I,J) ,J=1,8) 
101 FORMAT (A4,7F8.4) 
1 CONTINUE 


M=N-1 
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DO 12 I=1,M 
IF (WE(1+1,,2) Ge wrir.2))eo 10-12 
DO 13 J=1,8 
S(J)=WT (I+1,J) 
WT (I+1,J)=WT (1,7) 
WT (I,J) =S (J) 
CONTINUE 
M=M-1 
IF (M.GE.1)GO TO 3 
DO 40. I=1,N 
GQ=WT (I,2)/100.0 
SWT (I) =(GQ**2)*100.0 
DO 20.=1,N 
WRITE (6,201) SWT (I), (WT(I,J) ,J=1,8) 
FORMAT (/10X,F8.4,5X,A4,7(5X,F8.4) ) 
CONTINUE 
STOP 


END 
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APPENDIX 3 
RANK ORDER MATRIX OF GRAIN SIZE PERCENTAGES 
Sediment Grain Size Classes Grouped into the 


Type Wentworth Grain Size Scale 
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L 99-51 0.320 0.16 Oau2 0.0 0.0 


L 99.52 0.46 0.03 0.0 0.0 0.0 
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L 99:. 54 05.33 Oe) 02035 0.0 0.0 
L 99.54 0.42 0.04 0.0 O20 0.0 
L 997.56 0.44 050 020 020 6.0 
L 99°. o.0 0835 0.06 0.0 0.0 0.0 
L 99.64 0236 0.0 0.0 0.0 0.0 
L 99.64 On29 0.06 0.0 0.0 0.0 
L 99.86 0.14 0.0 0.0 0.0 0.0 
L 99.86 0.14 0.0 0.0 0.0 0.0 
L 9288 On LZ 02.0 0.0 0.0 0.0 
L 99.86 O12 0.0 0.0 0.10 0.0 
L 99°. 39 Deel: 0.0 02:0 0.0 0.0 
L wpekake Ns 0.09 OO 0.0 0.0 0.0 
L 99s OL: 009 0.0 05.0 050 0.0 


The symbols in the column headed 'sediment type' 


indicate the following sediment types: 


L Lacustrine sediments 
LT Laurentide till 

Cr Cordilleran till 

A Ash samples 

FS Fluvial sands 


Af Ice wedge samples 
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APPENDIX 4 
PRINCIPAL COMPONENT ANALYSIS 


The data consisted of ratio numbers or percentages. 
Tanner (1958) points out that cumulative percentage grain 
Size curves, which often "zig-zag", are only approximations 
of Gaussian curves, and that the irregularities point to 
forms of independence of different classes or groups of 
classes of grain sizes. In addition, each curve is 
independent, although it is known that certain curves will 
show distinct similarities due to common genesis. In this 
study the unknowns were the interrelationships between grain 
size classes, whilst the objective was to test for 
(pre-determined) groups of curves. In this case no 
hypothesis concerning the number of components or 
combinations of classes characteristic of each group was 
known. Since the form of the data suited principal component 


analysis this technique was used. 


The account of principal component analysis is taken 


from Kingy Be?s77 PICs 
The solution to principal component analysis is 
essentially based on a geometric manipulation applied to 


matrix algebra. The object is to produce a new set of 


variables y., which are uncorrelated one with another, 
i 
from an original set of variables x,, which may well have 


been highly correlated. The solution involves an orthogonal 
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transformation of Xs to produce Yar such that the total 
variance associated with x is preserved in the components 

of Ys The geometry basically proves that the matrix [U] 

and its transpose (uj? when multiplied, equal the identity 
matrix, such that [x] = [U] + [y], and hence [y] = [uj rt lhe 
Lawley and Maxwell (1963), and Anderson (1958) show that 

this property of an orthogonal matrix may be manipulated to 
give a matrix containing the variance and covariance terms 


for the y's in the expression 


where © is the variance-covariance matrix. Since the y's are 
assumed uncorrelated the off-diagonal elements of the multi- 
plies expression are the covariance elements, and equal to zero. 
Hence, the diagonal will contain the variance of the y's which 
are usually referred to as A's. For example in the three 


variable case —- 


jo 0a 
ful” - £ + (ul =|0 A, 0 | = [Al 
0 0 dg | 


It can be proven that fee st ely os [vu]? in which 
case the general form of the matrix may be referred to 
instead of the geometric form. Since the variances on the 
diagonal are assumed to be unity in the transposed matrices, 


the solution to the above expression in terms of the variance 


(o) and the A's is 
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O117A = Oa 031 
O12 Oo2o-A FOr = @ = | 2-AI| 
Ger 023 O337A 
where I is the identity matrix, ana |z - AI] = 0 is the 


Characteristic equation of the matrix. The solution of the 
equation is a polynomial expression of an order equal to 
the number of variables (or components of the diagonal). 
The X's are therefore the latent roots or eigenvalues, 
which are obtained from the coefficients of the polynomial 
terms in the characteristic equation. Since the eigenvalues 
in effect represent the amount of variance explained by each 
component, the rank ordering of the eigenvalues in the 
eigenvector will produce a canonical component listing. By 
substituting the eigenvector into the expression for y, the 
y values or loadings may be computed. 

In the case of nonstandardized data each y value has 
to be scaled by the corresponding square root of the 


eigenvalue. In this case the solution for y becomes 


for an expression in terms of the variance-covariance matrix. 
By matrix manipulation and adoption of the convention of 

L 

1 ; 
Lawley and Maxwell (1963) of (W) = TU) -< EP, ehenett can 


be shown that the characteristic equation becomes 
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where R is the matrix of correlations, which relates to the 


model 


where [Z] is the vector of standardized original variables. 
The latter is most widely used because of the 

computational efficiency. It has the disadvantage of not 

having any means by which the significance level of the 


results can be readily tested. 


-eolds.imsv isaiprixo sips Aiasbasd..20 
eft 20 seunged boew yiebkw ’s hai 
joa Io spstasvbse th eas asd eae OF 


edd to Level sonsoitinpse edd 


‘an et ae . pashan fs x f r 
NG m3 ' ie xa iy iq : , 
| | " ee 1 ca ak 
Ce hte anita Met 


“iii sy 7 gin 
et a 
iv fbgdh-caghcaie ea 
adi tae b hi rade apthonine te ie mers rie 
danlres Ln HonoRa tes ee bane: 

het ie hgh wit, Cy Ape bs 


ane a aie a 


; 
aun 


ns ae rs ria ; inne J 
‘ ¥ 7 y¥ , . : 5 > j > ar 7 ‘ . ‘ : Die , 
f wy _ iu ' 7 oh te U) a ; fi hey f , , ie) pa aan ie tr } 7 


167 


SET Sd To oF vate an bs a L£8T Sd eer 6€ ss i an 
£9t to So ot v6 o EO a, S9L ~Su cv? = BS Bla a ae 
97T LO Se: 4 VIL “oe 87 °T L8 sa Sh"? 9ST ae i ll 
99 sa 6O% e9t “oe Se 69 sa he ae O6T V 09°? 
Sol ao Or’ Et. Sa Cec L6 Sa ess cyVt Lo Eo Y 
31) oe Ve € 4 sa Soe 80 Sa bo ¢ T?t tod 8o° 7 
6 LO La ¢ L6 Sa Ve ‘ Sa Ve.-s 6 LO ce 3 
vot ou 9S -e Sal 4S psy ae i Sa veo Ss OLE 34 ES 9 
s sa BO°? “ sa vec jaye “i Ses é9t Stl EVO 
ee Su a tet Sa oS 39 sa oh A 6€T wo Ns iy 48 
L6 sa OV € Pet “out i al di sa res ct Sa oh: sae 
¥ Sa OS *E OZE Sd al SR O8T Sd 68° S fer V ec el 
v c G ni 


*a°r odAy jueuoduog «qtr eddy Queuodwop ‘qr edAy 3zusuodwuop “qtr edAz, 3ueuOoduwo5 


SISATVNY LNHYNOdWOO 'TVdIONIYd WOdd GALVYHNED SONIGVOT YOLOVA AO XIULVW 


G XIGNdddV 


= 
f 


hy 
» 


se ES ape 


é 
me 


oP tome a 
E 
. 


Sat 


168 


COT 
68 
SvT 
OLT 
Slt 
L8 
9¢ 
6ST 
LUT 
CUT 
eval 
O97 
O8T 
69 
LS 
T6 
EZ 


6€ 


Ef 4orls) 
t3a0 
78°0 
68°0 
LG 0 
G6°0 
L6°0 
vO°T 
Cite tl 
Scr b 
VoL 
62° 
Sit 
L9AL 
Co* 1 
Cora 
S9*a 


B9*¢ 


69T 
78 
98 
Sv 
SdT 
97 
8TT 
G9 
TS 
v9 
is 
OTT 
LET 
€6t 
£61 
88 
Eé 


vITt 


86°0 
86°0 
66°0 
66°0 
66°0 
66°0 
66°0 
66°0 
66°0 
66°0 
66°0 
00* at 
QO*e 
00°T 
LO te 
LO +a 
E-* 


TO * cf 


Cct 


£StT 
ECT. 
oF 
OL 
& 
6€ 
6ZT 
OLT 
vOT 


78 


Sa 


LO 


LO 


VEE 
EvASG 
80°C 
EG*Z 
6S*<¢ 
Skee 
Lot? 
L8*¢ 
Loe Cc 
eG*e 
SU. 
GA 
pore 
Este 
Sar 
SBre 
S6°E 


Ere? 


Get 


m m 
So im 
a 

bar 

Sa -< 
ae] m 
“i 2) 
° * 
Oo 
fo A 
a @ 
ee 

f im 
oO cB 
fia) en 
. ’ 
© oS 
” yv 
"> bow 
— 

=) te 
tu sad 
Ne 
* ” 
ee 
at 
a” 
4» @ 


= = 
> 


GL J¥2 


0*S¥ 
0°89 


T°A0 


B2 TAO 


ra 


GL YAP 


o-aT 


T*se 


‘CL 


0° a2 


rsa 


a3 


PY sre 


9°33 


“43 ke 08 
3" 2 PL 3 


rie 
ssi Gl Tre 


169 


penutjzUuoostp 
uoTRZeHbTASSAUL 


OST Lo O20 


vv 
8S 
cL 
OT 
SS 
8T 
6€ 
€9 
SET 
LST 
es 
G8 
bye 
SOT 
VEL 
Stt 
89 


c8T 


76 °O 
€6°0 
€6°0 
€6°0O 
€6°0O 
€6°0 
v6°0 
96°0 
96°0 
96°06 
96° 
96°0 
96°0 
£6*¢@ 
£6°0 
L6°0 
86°0 


86°0 


EST 


SES 


VST 


O9T 


BET 


8LT 


62 


60 


64 


di 


ice 
ECO 
dean0 
€7°0 
Eaa0 
85-0 
v.90 
GE +0 
C820 
98°0 
96*9 
£0°T 
S$O°D 
et*s 
eT a 
OG *E 
Be°t 


6V°T 


L8 


eer 


8cT 


O8T 


C6T 


Tél 


OVT 


€vt 


6vT 


VOUT 


88T 


68T 


Est 


L6 


87T 


Sa 


LO 


T6°0 
S6*0 
€6°0 
66°0 
60°T 
So" oD 
ttt 
vv'T 
OS°T 
¢S*t 
ES *t 
69° 
O4£*E 
Gs *£ 
96° a 
6L°<¢ 
GS" a 


CG *S 


ad 


0° aS 


T23 


r 


r 28 


0°33 
0°33 


32 
38 


roe 


AS 


ra 


TO 
33 


0°33 


0°33 


- 


eR ERS EGRSS 
che Om geen 


170 


CY 
SCE 
y 
cg 
vS 
DLT 
8L 
mS 
OS 
SET 
a 
V8T 
OOT 
£9 
C8 
9G 
Fe) 


98T 


L'T 


L'T 


ca @ 
CLO 
jy ls) 
ae O 
CL 0 
(ay 3) 
8L°0 
08°0 
LeU 
v8°O 
S8°0 
98°0 
88°0 
06°0 
06°0 
T6°0 
joel, 


c6°0 


poenutjUuoostp 
uoTRIeHTASSAUL 


e9L At 7T°O 


[ene Shak Sits0 
a) oe I al Lbs 0 
€6 Sd 0c°*0 


LY LT be-20 


6LT 


O9T 


po a 
OEE 
T6T 
ce 
VoL 
St 
06 


Scr 


sa 


LT 


Sa 


Sa 


LT 


LT 


sad 


LO 


Sa 


LT 


LO 


LO 


LO 


LT 


Sa 


LT 


LT 


LO 


Ma 
Gc 0 
VEoU 
87°0 
Sen) 
Soe 
Veo 
ve°o 
Lot 
Srau 
ape eat) 
ae Ea 
Oana 
vL*O 
St bi 0) 
90 
Cow 


ts 0 


ye 


aT Sale 


a a 


ye 


gol 
sel 


gh 


ria 


02 
ee 
sy 


o 
ee] 


Le 


a 


+ 
A 


Ta 


= 


* 


a < = 
Tl 


ul 


ba 


r 


38-0 
- : 28.0 - 
7 $8.0 2 


a4 
08.0 
8y.9 
Te.0 


vv.0 


vt -0 


«88 -0 i 


Ce ne 


oe 


ee 
ear ee ae 
ef 


an 


2 
_-2 pe s5 
~<a oa 


ee een Gh 


Lat 


6S 

881 
8ST 
DL 

O6T 
OcT 
COL 
Bek 


eer 


Al A 


LO 


seo 
cv °O0 
evo 
cv 0 
Lv°0 
87°0 
OG 0. 
SSG 
Sce0 
SO) 
0:9°>0 
09-0 
v.90 
S90 
99" 0 
990 
OL°0 


LL’ 0 


oe 
09 
EGU 
est 
vE 
82 
sé 
SOT 
L8T 
€6 
SLL 
76 
8ST 
O€ 
VC 
9% 
COT 


nd 


20 
782 
28 

Aq 
T30 


ehF eee 
i s + Ms te e a 


Stites 


® 2 ¢ 


Cpe a 
ee A, saa 


ca 330 
PL. Tos 


A3 
ra 


, 
r 
r 
r 
¥ 


0°32 
O° #5 

"€3 
Oo" €5 
0° ¢3 
o" #38 
0°26 
0°23 
0" 22 


- 7 
li Tin 


, or) ea 
jihaeers eet ce 


a as 


sagezega 
piscanere 


rh 


Shedd 3 = | ae 


172 


60T 


LT 


LT 


LT 


L'T 


c0°O 
€0°0 
vo°o 
L0°*0 
L050 
OT°O 
PE S6 
20 
Le 
LES 
SiO 
8T°0 
6T°O 
ca 0 
qo 0 
CE"0 
Ero 


ve°o 


Sst 
Eel 
LOT 
EOT 
TST 
Ge 
LY 
8E 
Lit 
02 
VT 
gt 
90T 
68 
OL 
66 
&¢ 


Se 


LT 


L'T 


LT 


LT 


L'T 


Sa 


so: U= 
66 0= 
£6. 0° 
66° OQ" 
v6 0= 
v6 0= 
co US 
Eo 0" 
Ge 0- 
ts‘ O= 
08 O= 
68° 0= 
0s 0= 
Sa 0= 
va G= 
ec’ = 
Eo = 


vo 6" 


«3 
2 
38 
3 
33 
; 


z78 


penutjuoostp 


uOoTIeHTASSAUT 


6L 
9S 1 


Sit 


86 


LT 


sa 


os 'L= 
9e° i= 
UE i= 
oe” i= 
3G L= 
eée° l= 
s¢° = 
eo t= 
ee l= 
Oc T= 
Sit ill Sag 
30° {= 
&0° = 
60" y= 
50 t= 
Zo" = 
vO" t= 


vO L= 


1 . ; 
‘geal 
4 <2 ; | | | 
j @ “ : | : ; | 
p by 4 i uae ee i 
eo , aes aid ‘ ~~ ny, fi J ‘a sie ry 
| ie | P Ay) ih ona . 7 a. ‘ie ai | i 
pS pt | ; | 
J LY ted ; | ; 
f oD ; ; | 
a i an ' 7 t | 
| a ; : 7 1 , 6 : 
ah : | | 
| | | TAL 
( 
= 


ty 
ree 
a) 


46 ‘ 
igh = 


} ; 7 


+e ) 


a 


as 


te ae 


iy 


174 


£9 


v8T 


CrlaL 


EST 


VST 


89*t= 
po<T= 
EO LT = 
CO ed = 
L2-T= 
PO ' T= 
£S*t> 
9S*4= 
96+ t= 
9G* t= 
gee 
cS* t= 
67-i= 
‘pet= 
So* t= 
SD fae 
iv * ts 


BS" ts 


173 


96 
r9 
L9 
(Ae) 
VOT 
98T 
vL 
OOT 
SET 
vt 
Cv 
ss 
T8 
Tv 
89T 
CTT 
8L 


eS 


oo" i= 
56° [= 
68 T- 
Ge c= 
uo. T= 
Se t= 
30° t= 
ea 
Eo t= 
V6 = 
eo [= 
sé°T- 
te T= 
96° T= 
va t= 
6° T= 
69° T= 


69° T= 


176 


OTT 
eo 
97 
jy 
S8 
OT 
OS 
VEL 
78 
ee 
Co 
8T 
cs 
vO 
IST 
eS 
8S 


9S 


EG" r= 
70° :E= 
ol he 
a ea 
oes des 
C1 a aa 
36> b> 
SO t= 
So t= 
VO ne 
ro c= 
Lot 
VON 
Uae no 
Coe 
CO: f= 
Cot 


tour 


* ‘ : 7 | i. ho, vuetd 
an btieiel dice 3 eh te o ne a 


L779 


T6T 
VEL 
v6 
COT 
88 
8TT 
69T 
cs 
St 
Le 
se) 
v9 
SET 
89 
SoG 
TT 
98 


eit 


66-1.= 
So-1= 
66 t= 
6Gu,— 
66. 7= 
662 -.= 
667 = 
6GmeE= 
eo ie 
86° t= 
SOs = 
86r- T= 
S61 - 
S6r= 
She t= 
UCaar— 
£6 T= 


G6 t= 


178 


seTdues shpam sof I 
sotdues usv W 

so—Tdues pues TeTAnTg Sua 
sorTtdues [[T} uexreT[Ttpsz09 To) 
sotdwes [TT eptjuezney GT 
sotdues sutizqzsnoeq T 


:sodAj jueuTpes HuTMOT TOS 
24} JUeSeadez uUMTOO sTy} UT USsATH sepod syuL ‘*pTeTz eyy Lin 
opeuw uot zeooTTe ed&j WuosulTpes TeuThtz70 |yA sejzeoTtput rwedfz,, 


ST3T2 SUR YRTM UUNTOO |yW ‘/xTAAeu Hutobezo0Z ayy ut 


Se 


oiab 


ok 


i 


667 7E= 


66 “p= 


asiqmse [fis ebitnex 


selamse ILis nexellibiod 


/ aslqmse buse Isivult 
seni eslqmse eA 
eslqusa sphsw sol 


APPENDIX 6 


PROCEDURE FOR TREND SURFACE ANALYSIS 


DATA COLLECTION 


The assumption made concerning the data to be used 


in the analysis was that summit peaks in the Elbow River 


basin are remnants of an earlier cycle of erosion. The data 


therefore represent the highest peaks within each positively 


Closed contour in the basin. Two criteria were, however, 


applied before the final point value was utilized; 


iy) 


2) 


If two points were closer than one inch, the highest 
value was taken. (The data was collected on 

1:50,000 topographic maps). 

The summit peaks considered ‘ere selected so as to 
be most representative of a summit surface if they 


occurred on valley divides. 
The Analysis 


"A special case of the multiple regression problem 
- » - is trend surface analysis. For a map pattern 
of a variable Z, which is measured at a number of 
discrete location points, the analysis involves the 
regressing of this dependent variable on two 
independent variables U and V, which are orthogonal 
coordinates for the map area." 


King, Ue J. “Gl969).)) jp. 1S7 0 


As Krumbein and Graybill (1965) point out, the Z 


values are divided into two parts by the analysis. The first 
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part is associated with the "large scale systemmatic changes 
that extend from one map edge to the other. . .". This is 
the trend surface, or surface of the regression. The second 
part is associated with ". . . small-scale apparently 
nonsystemmatic fluctuations. . .". This refers to the 
deviations of each point from the surface, and concerns 
the residual values. Board (1965) refers to part one and 
part two as regional and local effects respectively, whilst 
Chorley and Haggett (1965) refer to these two parts as signal 
and noise. It is these two components that are used to 
interpret the data. 

The surface is constructed such that the sum of the 
squared residual values is as small as possible using the 


basic equation of the model 


Zi5 = TAU ae) ~ See (King) be Js) 1960) 
where Z is the observed value (elevation) at location 
(g,V). The trend is TAO and oes is the residual. 

To calculate the trend parameters a polynomial 


expression is used, and in this#@case calculated to the sixth 


order which involves 28 terms. The form of the polynominal is 


Z Z 
. . . + eee#ee + . . 
U,+a V.+0 U; +05,U;V +O poV e 
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The computations involved when the original data is 
taken on a regular grid pattern is simplified by the use of 
orthogonal tables. In this application the points were 


irregularly spaced and the computations more complicated. 
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Whitten (1963) described methods which allow easier 
computation, especially for irregularly spaced control 
points. The computation involves a summation of each of the 
terms of the polynomial. This procedure is now commonly 
used in the programs available for trend surface analysis. 

The program used (O'Leary, Lippert, and Spitz, 1966) 
in this study, utilized the same technique for computation 
of the trend as that outlined by Whitten (1963), for 
irregularly spaced control points (Krumbein, 1959). The 
following are the computations used for the error measures 
in the computer program. 


Total variation (V) 


N 
G eaproeb(gne-iz) ¢ 
=) 
N is the number of observations. 
N 
i=1 


al: 
ae is the calculated value. 


Variation explained by the surface (E) 


pe oae so 


Coefficient of determination (T) 


T = E/V 
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Coefficient of correlation (C) 


Standard deviation (d) 
L 
by a Syn) 


The goodness of fit of the surface was tested using 
analysis of variance. The method is clearly described by 


King; C. A. M. (1967), and uses an 'F’ test. 


Between group variance 


F value = Within group variance 


where Q = S/N-(df-1) and R = E/(df-1) and df is the 
number of expressions in the order of polynomial used to 


derive the surface concerned. 
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